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Acesta este prima parte a cursului predat studentilor sectiei de Automatica a Facultatii de 
Automatica si Calculatoare, Universitatea Politehnica din Bucuresti, in anii 2002-2003 
si 2003-2004. Textul reprezinta in fapt note de curs (folii) ce au fost videoproiectate 
studentilor. Orice sugestii si corectii sunt binevenite | 


Corespunzator impartirii in doua semestre, cursul cuprinde doua parti principale: Teoria 
Sistemelor | si Il. Detaliem in continuare structura cursului TS | pe capitole si sectiuni: 


1 INTRODUCERE 


1 Conceptul de Sistem. Exemple 
2 Abordari Fundamentale 
3 Sisteme de Reglare Automata 


2 SEMNALE SI SISTEME 


1 Semnale 
Definitie. Exemple 
Clase de Semnale 
Operatii cu Semnale 
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Semnale cu Impulsuri (Semnale Singulare) 
2 Transformari Integrale 
3 Sisteme. Proprietati Fundamentale 
4 Sisteme de Convolutie 
5 Raspuns in Freceventa. Functie de Transfer 
6 Sisteme de Convolutie cu Functia de Transfer Rationala 


3 SISTEME CONTINUE SISO 


1 Raspunsul sistemelor SISO 
2 Stabilitate 
3 Regim Permanent si Tranzitoriu 
Regimuri de Functionare Fundamentale 
Raspunsul Sistemelor la Semnale de Intrare Standard 
Performante de Regim Permanent si Tranzitoriu 
Raspunsul Sistemelor Elementare la Intrari Standard 
4 Reprezentarea in Frecventa a Sistemelor 
Diagrame Bode 
Diagrama Nyquist (Hodograful) 
Criterii Frecventiale de Stabilitate 
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4 CONCEPTE FUNDAMENTALE ALE BUCLEI DE REACTIE (FEEDBACK) 


1 Ce Dorim ? 

2 Bucla de Reactie 

3 Stabilitatea Interna a Buclei de Reactie 
4 Stabilizare 

5 Performante Asimptotice 

6 Solutia Problemei Reglarii 

7 Exemple 
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CAPITOLUL 1: INTRODUCERE 
Motivaţie. Automatica în cadrul științelor inginerești. 


Ingineria clasică (pâna la al doilea război mondial): 

- dispozitive/instalaţii de transformare a energiei (mașina cu aburi, motorul electric); 
caracteristici principale: randament, debit. 

- comanda/acţionarea unor astfel de instalaţii: operator uman. Informațiile nece- 
sare acționării parvin prin intermediul simțurilor umane, în urma observaţiilor privind 
rezultatele acţiunilor anterioare. 


Ingineria moderna: dezvoltările tehnologice au creat dispozitive a căror acţionare directă 
de către om devine imposibilă datorită limitărilor fiziologice; 

Exemplu: Construirea unui sistem automat de tragere antiaeriană. 

Motivaţie: la sfarșitul celui de-al doilea război mondial, viteza țintelor a devenit 
comparabilă cu cea a proiectilelor, astfel că tirul antiaerian nu mai putea fi comandat cu 
succes doar de către operatorul uman. 


Figura 1: Schema simplificată a ghidajului unei rachete antiaeriene 
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Elemente necesare pentru funcționarea unui astfel de sistem: 
- informaţia asupra situaţiei momentane 
- extragerea datelor utilizabile pentru comandă 


- manipularea cestor date în vederea elaborării comenzii, având un scop bine determinat: 
doborârea țintei 


Problema se reduce la construirea anumitor dispozitive de comunicații și comandă, 
pentru a căror proiectare sunt necesare elemente de 


- teoria comunicaţiilor 
- tehnică de calcul 
- teoria sistemelor automate 


Norbert Wiener, 1948: “Cybernetics, or Control and Communication in the Man and 
the Machine” 
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Instalaţiile din ingineria clasică: esenţial este randamentul - aspect de tip energetic. 


Dispozitive de comunicaţie și comandă: esenţială este acuratețea cu care este transmisă 
și prelucrată informaţia. 


In proiectarea dispozitivelor de comunicaţie și comandă se face abstracţie de natura 
(mecanică, electrică, chimică, etc.) procesului /instalaţiei/ sistemului tehnologic. 


Exemplu: un receptor radio primește o cantitate foarte mică din energia emiţătorului, 
dar important este ca semnalul recepționat să fie “curat” (fără zgomot sau distorsiuni). 


Teorie generală a SISTEMELOR. 


Noţiunea de SISTEM: conceptul fundamental. 
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1. Conceptul de Sistem. Exemple 


DEX: “Ansamblu de elemente (principii, reguli) dependente între ele și formând un 
întreg organizat, care pune ordine într-un domeniu de gândire teoretică”, reglementează 
clasificarea materialului într-un domeniu de știinte ale naturii“*, sau face o activitate 
practică să funcţioneze conform scopului urmărit***." 


* Sisteme social-politice; sisteme filozofice. 


xX 


Sisteme fizice; sisteme biologice. 


*** Sistem informatic; sistem energetic; instalație tehnologică, proces. 
Sistem (in inginerie): “Ansamblu bine organizat de părți interdependente, capabil să 
răspundă, sub acțiunea a diverși stimuli, unui anumit scop, cu anumite performanțe” 


Exemple: un automobil, un sistem de operare, un sistem informatic. 


Capitolul 1 - Introducere. Conceptul de Sistem. Exemple 


Teoria (matematică a) sistemelor (automate) 
Noţiunea de sistem: un concept matematic precis. 
Definiţia 1 (Kalman, 1969). Concepem un sistem ca o structură în care ceva (materie, 


energie, informaţie) poate fi introdus la un anumit moment dat și din care rezultă spre 
exterior ceva la un alt moment de timp. 


Figura 2: Model de tip “BLACK-BOX" 


e u si y (“ceva”) sunt semnale. 
e Problema: ce caracteristici are un astfel de model ? 


e Sistem: “operator”. 
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Exemple 


1) Amplificatorul operaţional. 


Figura 3: Amplificator operaţional 


vo = A(v4 —v_) E Av. 

Transformă un semnal în alt semnal: v_— vo. 

Q: lărgimea de bandă, răspuns la intrare sinusoidală? 
2) Ecuația diferenţiala t = f(x), z(0) = zo 
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Conceptul de Sistem. Exemple 


Asociază o soluție unei condiţii inițiale date: £o — ¢ġ(t; £o). 


Q: creștere exponențială, periodicitate ? 
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Exemple 
3) Program de calculator: transformă un șir de caractere în alt șir de caractere. 


De exemplu, auto — otua. 


Q: regula ? 

4) 
f tii tim u 
S : ©- ys Iu 
" m tpm Um 


Unui vector din R’” i se asociază un alt vector din RP. 
5) 
T 
„= | true) ar 
0 
Unui semnal cu suport în |0, T] i se asociază un alt semnal cu suport în [0, T]. 


Q: Ce au in comun sistemele definite la 4) și 5) ? 
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Divizorul de tensiune în curent alternativ 


vi R = Vo 


le, O 


Figura 4: Divizorul de tensiune în curent alternativ 


= | dt FLS 
0 T | v = Ri 
Comportament Intrare/leșire (1/0) Comportament Intern 


Proprietate importantă : Liniaritatea 


R și L sunt elemente liniare de circuit. 
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2. Abordari Fundamentale 


SISTEM: “operator” care transformă intrarea în iesire. 


S-au impus două puncte de vedere: 


a) definirea noţiunii prin caracterizarea comportamentului intrare-ieșire = maniera 
operatorială de abordare. 


Analiză funcţională 


b) elementul primar îl constituie comportamentul intern = maniera newtoniană de 
abordare. 


Ecuaţii diferenţiale 
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Modelare matematică 


Modelare matematică: procesul de scriere a unor ecuații care stabilesc dependenţa dintre 
diverse mărimi care acţioneaza în sistem sau asupra sistemului. 


a) pe baza datelor experimentale (intrare-ieșire) din sistem: Identificarea Sistemelor. 


b) pe baza ecuațiilor fizico-matematice: abordarea newtoniană. 


Exemplu: dependenţa tensiune-curent a unui rezistor. 


a) î = 1.02 1.98 3.03 4.01 5.01 mA 
v = 1.01 2.00 2.99 4.00 5.02 V 
—> v (aproximativ) proporțional cu i: v = Ri; rezultă R = 1k0. 


b) Legea lui Ohm 
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In elaborarea modelului este necesar un compromis între complexitatea modelului și 
descrierea adecvată a fenomenului / procesului respectiv. 


Model satisfăcător: diferenţele între rezultatele obținute prin calcul și cele obţinute 
experimental trebuie să fie mai mici decât anumite toleranțe impuse. 


SCOP: Proiectarea Sistemelor de Reglare Automată (SRA). 


Capitolul 1 - Introducere. Abordari Fundamentale 


3. Sisteme de Reglare Automata 


Reglare automată: procesul de a impune ca anumite variabile specificate ale unui sistem 
să urmeze anumite evoluții impuse, în prezența diferitelor perturbații (precum și a 
incertitudinilor de modelare). 


Influențarea evoluţiei unui sistem fără intervenţia umană. 
Qp 
Qout 


Tref 


Termostat ~~ Valva gaz Cazan EO Proces £ 


Qin 


Figura 5: Schema bloc de reglare a temperaturii într-o casă 
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Diagrama unui Sistem de Reglare Automată 


| 
Perturbatie | 


Referinta Eroare Iesire 


> — | Regulator Element | y Proces a 
B de executie 


Iesire masurata 


Senzor 4 


A 


Zgomot 


Figura 6: Diagrama bloc funcţională a unui SRA 


Observaţii: 
- ansamblu de sisteme interconectate 
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- existenţa reacției inverse | 


Reglarea cu reacţie inversă folosește măsurători ale variabilei reglate pentru a influența 
mărimile de intrare ale sistemului reglat, astfel încât variabila reglată să urmărească o 
anumită evoluție impusă . 


Conexiunea inversă permite evaluarea permanentă a gradului de realizare a obiectului 
propus: metodă extrem de eficientă (vezi exemplul următor). 


Pe parcursul expunerii, vom utiliza frecvent schema bloc simplifcată din figura 7. Această 
schemă prezintă configuraţia standard a unui sistem în reacție inversă. (In practica, se 
utilizează diverse scheme de reglare automată, care pot fi compuse din mai multe bucle 
de reglare /regulatoare, specifice fiecărei aplicaţii în parte). 


Figura 7: Schemă bloc simplificată a unui sistem în reacţie inversă 
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Obiectivul principal al cursului constă in descrierea unor metode generale de sinteză a 
compensatoarelor (regulatoarelor). Pentru a putea înțelege și aplica aceste tehnici de 
sinteză, este necesar mai întâ să studiem anumite proprietăți și caracteristici sistemice; 
în consecință, prima parte a cursului este dedicată metodelor de analiză a sistemelor. 


Capitolul 1 - Introducere. Sisteme de Reglare Automata 


Reglarea vitezei unui automobil 


O primă analiză (calitativă) a avantajelor reglării cu reacție inversă. 
Se pune în evidenţă rolul fundamental al reacției inverse în sistemele de reglare automată. 
Uref = 55km/h - drum orizontal. 


Factori de influenţă: 
- accelerația: unghiul clapetei de accelerație (măsurat în grade °) 
- panta drumului: gradul de înclinare (măsurat in %) 


lpoteze: 


a) avem regim staționar (nu există dinamică) - noţiunea va fi introdusă în capitolul 3. 


b) sistemul este liniar: 
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- la un unghi u de 1° al clapetei de acceleraţie, viteza crește cu 10km/h. 
- la o înclinare a drumului w de 1%, viteza scade cu 5km/h. 


v = 10u — 5w 
c) măsurătorile sunt precise (eroare de până la +1km/h) 


Utilizăm 2 scheme de reglare: în buclă deschisă și, respectiv, în buclă închisă. 


|. Reglare în buclă deschisă 


0.5 


| 
| 
u -= 
Lo 0.1 rA E 10 (IRAN 


Figura 8: Schemă de reglare în buclă deschisă a vitezei unui automobil 
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Notăm viteza de referință vref cu r. In acest caz, u = 0.1r și viteza automobilului este 
dată de j 
v = 10u — 5w = 10(u — 0.5w) = LOU) — 5w = r — 5w. 
Performanțe: 
Eroarea absolută: € := r — v = 5w 
Eroarea relativă: e, = e/r = 5w/r. 


w = 0: v = 55km/h, e = 0km/h, e, = 0%. 
w = 1: v = 50km/h, e = 5km/h, e, ~ 9.1%. Performanţă destul de slabă. 
w = 10: Erori f. mari! (exercițiu) 


Dacă “amplificarea” sistemului se modifică sau am estimat-o de la început imprecis (la 
un unghi u de 1° al clapetei de acceleraţie, viteza crește cu 9km/h în loc de 10km/h) 
cum se schimbă performanţa sistemului de reglare proiectat ? 


v = 9u — 4.5w = 0.9r — 4.5w, e =0.lr+4.5w, e, = 0.1 + 4.5w/r. 


Chiar și în absența perturbaţiei (w = 0) eroarea este semnificativă: e = 5.5Bkm/h, 
e, = 10%. Situaţia se înrăutățește dacă w 70. 
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Concluzii: 


- erori relative mari (chiar și la variații mici ale perturbaţiei) 


- erorile în modelarea sistemului (10%) se transmit integral în eroarea relativă (10%) 
(nu se atenuează deloc, nu avem robustețe !) 


- reglarea în buclă deschisă nu este o soluție “buna” 


Cum corectăm neajunsurile observate ? 


||. Reglare în buclă închisă 
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ii g u - v 
+ 
== mic ÎI ai, K=100 A > 10 2 


Figura 9: Schemă de reglare în buclă deschisă a vitezei unui automobil 


Analizând schema de mai sus (figura 9) rezultă că viteza automobilului este 
v = 10(u — 0.5w) = 10(100(r — v) — 0.5w) = 1000r — 1000v — 5w, 
de und 
S unas 1000 5 
v = —r — ZU 
1001 1001 


Capitolul 1 - Introducere. Sisteme de Reglare Automata 


Performanţe: 

Eroarea absolută: e := r — v = 0.001r + 0.005w. 
Eroarea relativă: e, = e/r = 0.001 + 0.005w/r. 
Ambele erori sunt “mici” ! Performanţă bună. 


w = 0: v = 54.945km/h, e = 0.055km/h, e, = 0.1%. 

w = 1: v = 54.940km/h, e = 0.060km/h, e,. ~ 0.11%. Creștere nesemnificativă în 
raport cu situația w = Q. 

w = 10: v = 54.895km/h, e = 0.105km/h, e, ~ 0.2%. La o creștere de 10 ori a 
perturbaţiei, eroarea relativă s-a dublat. 


Dacă “amplificarea” sistemului este 9km/h (în loc de 10km/h) ce se întâmplă cu 
performanţele sistemului de reglare în buclă închisă ? 


_ 900, 4.5 l 4.5 
i e = 4 = —r FU = lir + 0.0049 
= 0 901 = 0.9988r — 0.0049w, e = 901” rF LU > COULIE w, 


ep = 0.0011 + 0.0049w/r. 


Nu se schimbă aproape nimic! Erorile absolută și relativă sunt practic aceleași, deși am 
avut o eroare relativă de 10% în “amplificarea” sistemului. 
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Exercițiu: Analizaţi modificările apărute în determinarea erorilor, în cazurile w = 1 și 
w = 10. 


Explicație: Amplificare mare și reacție inversă! 


Concluzii: 


- în cazul w = 0 nu avem reglare perfectă, adică v Æ r; eroarea (de regim staționar) 
nu este zero, dar este foarte mică. 


- dacă amplificarea în buclă crește, atunci: 


- creşte precizia în regim staționar (scade eroarea de regim staționar) 

- scade senzitivitatea (relativă) a mărimii reglate în raport cu perturbația 

- scade senzitivitatea (relativă) a mărimii reglate în raport cu imperfecțiunile mod- 
elului. 


- reglarea în buclă închisă este eficientă. 
Amplificare mare ——> performanțe mai bune. 


Capitolul 1 - Introducere. Sisteme de Reglare Automata 


Q: Cât de mult poate crește amplificarea ? 


Nu poate crește oricât, apare fenomenul de instabilitate în buclă. (Nu putem ține 
accelerația “la blană”). 


Stabilitatea: cerința principală de proiectare 
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Proiectarea SRA 


1. Stabilirea obiectivelor reglării. 

2. Identificarea mărimilor ce trebuie reglate. 

3. Scrierea specificaţiilor pt. mărimile reglate. 
4. Stabilirea configurației de reglare. 


5. Obţinerea unui model (pt. fiecare element al buclei: proces, senzor, element de 
acţionare, etc.) 


6. Alegerea unui regulator și a parametrilor acestuia. 
7. Optimizarea parametrilor și analiza performanţelor. 


8. Dacă performanţe = specificaţii atunci proiectarea este încheiată. In caz contrar, se 
reia “algoritmul” de la pasul 4. 
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Obiectivele proiectării 


: specificaţii (ce trebuie realizat). 


Caracteristici: 

- complexitate 

- compromisuri 

- FISC 

- nu iese întotdeauna ce și-a propus proiectantul 


Analiză și sinteză ... 
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Compromisul de proiectare 


1) STABILITATE: cerința principală de proiectare 
2) PERFORMANȚA: 


- eroare de reglare (în regim staționar) mai mică decât o valoare specificată 


- dinamica răspunsului =—> parametri specificaţi 


3) ROBUSTEȚE: Stabilitatea și performanţele trebuie menținute în condiții de model 
cu incertitudini. 
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CAPITOLUL 2: SEMNALE SI SISTEME 


1. Semnale 
1.1 Definiţie. Exemple. 


Un semnal este o funcție de timp. 


Mărimi fizice variabile în timp: forța F care acţionează asupra unui punct material, 
tensiunea vo la ieșirea unui AO, presiunea p a unui fluid. 


Notaţie. 

e F, vo, p sau F(-), vol), p(-) se referă la semnal sau funcție; 

e F(t), vo(10.33), p(t — 1) desemnează valoarea semnaleleor la momentele t, 10.33, 
t= 

Definiția 2. Se numește semnal continual o funcție f : T — A, unde A este o mulțime 
dată numită imaginea (sau mulțimea de valori a) semnalului iar 7 este axa (sau domeniul 
de definiție al) semnalului. 

Dacă 7 C R (mulțime “continuă” ), atunci u este un semnal continual; în cazul în care 


wv] 


T C Z (mulțime “discretă” ) atunci u este un semnal discret. 
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Semnale eșantionate 


Semnale eșantionate: falk) = f(to + kT) (mai multe detalii la Sisteme Discrete). 
T se numește perioada (sau pasul) de eșantionare. 


fik) | fo 


f(k)=f(0+kh) 


Figura 10: Semnal eșantionat 
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1) Cursul leu-dolar. Axa semnalului: discretă; imaginea: R.. 
2) O secvenţă semi-infinită de biți: 0111001.... Axa semnalului: Z4}; imaginea: {0,1}. 
3) Tensiunea de ieșire a unui AO. Axa semnalului: R4}; imaginea: R. 

) 


4) Nivelul apelor Dunării: semnal eșantionat. 
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Notaţie. Clasificări. 


e Notăm cu SA, SU mulțimea semnalelor continuale, respectiv discrete care iau valori în 
mulțimea A. 


e |n mod uzual vom lucra cu semnale reale sau complexe: imaginea semnalelor va fi 
A = R sau A = C. Mulţimi de semnale întâlnite frecvent: Sa, Sc, Să, Sa; 


e Semnalele pot avea 
a) axă finită: 7 = (a,b) - interval finit; Z = {n,n + 1,...,n +1} - mulțime finită; 


b) axă semi-infinită: 7 = R4 (t > 0), T = R- (t < 0) sau 7 = Z, (k > 0), 
Teha oD 


c) axă infinită: 7 = R, T = Z. 
e Semnale scalare (A = R sau A = C), semnale vectoriale (A = R” sau A = 0). 
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Semnale standard 


J4 20 
10 ż<0 
Semnal scalar real, continuu pe porțiuni. De tip “curent continuu”. 


a) Treaptă unitară: 1(t) 


> 
b) Treaptă unitară discretă: 1(k) = | | k>0 


t t20 


c) Rampă: ramp(t) = | 0 +<0 


Semnal scalar real, continuu pe porțiuni. Observaţie: ramp(t) = t1(t). 


Funcție de tip “polinomial”. 


In mod similar se definește și semnalul rampă discret. 
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1(t) 


ramp(t) 


Figura 11: Semnal treaptă și rampă 
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d) Impuls discret:  6(k) = | 0 kk F 0 


Figura 12: Impuls discret 


< < 
e) Impuls dreptunghiular: rect(t) 2 | l o b 
f) Impuls triunghiular: trian(t) = | 1 — jti h t<] 
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rect(t) 


-1/2 


trian(t) 


1/2 
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g) Semnal de tip armonic: 
u(t) = Acos(ut + $) 
A - amplitudinea 
w - pulsaţie; ww = 2r f = 27 /T unde f € R, este frecvența semnalului iar T € R} este 
perioada acestuia. 
p - faza (sau defazajul) 


Reprezentarea complexă a semnalelor armonice (a € C): 


u(t) = ae = Ae'?e'“i = Acos(ut + d) + j Asin(ut + e) 


Semnale “reale”. 

Radio (AM,FM) 

Satelit, cablu TV 

Video (Pal/Secam, NTSC) 
Telefonie (mobila, fixă) 


Aceste semnale nu sunt definite de formule matematice, ci de anumite caracteristici: 
frecvență amplitudine, factor de umplere, etc. 
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1.2 Clase de semnale 


Am menţionat deja că vom lucra frecvent cu semnale din Sp, S, Sc, Sâ. Aceste 
mulțimi de semnale sunt foarte bogate; vom considera (din raţiuni tehnice dar și 
utilitare) submulțimi ale acestora. In general, semnalele continuale sunt funcţii continue 
pe porţiuni și/sau local integrabile. 


Observaţia 3. Pe Sk și SË (Sc și S£) se poate introduce o structură de spațiu vectorial 
peste R (peste C). 


Clase de semnale: cu acţiune finită , de energie finită , mărginite. 


a) Semnale cu acţiune finită. 


+00 
L'(R) = {u € SR : u masurabila, J lu(7)|ldr < œ} - semnale absolut integrabile 
pe R. Bi 


+00 
2AE ` jz(k)| < œ} - semnale absolut sumabile pe Z. 
k=—00 
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Terminologia. Să considerăm un punct material de masă m asupra căruia acţionează o 


t2 
forță F, pe un interval dat de timp, finit sau infinit; avem J F(r)dr = m(v(t2)—v(t1)), 
ti i 
1 2 
unde v este viteza punctului material. Dacă v(t1) = 0, rezultă v(t2) = | Fa (cai le cai 
m Ji 


“Acţiunea” forței are ca efect modificarea vitezei punctului material. 


b) Semnale de energie finită. 


+00 
L? (R) = {u € Sg : u masurabila, J lu(r)| dr < oo - semnale de pătrat integrabil 


— OO 
pe R. 


+00 
(2) = (zei : > lz(k&)|2 < oo? - semnale de pătrat sumabil pe Z. 


k=—00 


Terminologia. Puterea unei rezistențe R parcursă de un curent de intensitate i(t) este 
P(t) = i2(6)R. Energia disipată pe un interval dat de timp, finit sau infinit, este 


f Poar f Poar 


tı ti 
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c) Semnale (esenţial) mărginite. 
L*(R) = {u € Sg : u masurabila, essup|u(t)| < oo! - semnale esențial mărginite pe 
R (mărginite aproape peste tot). 


Vom considera in mod uzual semnale mărginite: 


u E€ SR este mărginit dacă există M > 0 astfel incât |u(t)| < M (< oo), pentru orice 
teR. 


Orice semnal mărginit este și esenţial mărginit. Există semnale esenţial mărginite care 
nu sunt mărginite ? 


+00 
In cazul discret, 1° (Z) = {x € SI : > sup z(k)” < oo) - semnale mărginite pe Z. 
k=—00 
Observaţia 4. Vom considera în mod obișnuit semnale având axa de timp semi-infinită 


I*(Z,), LI(R4). Axa de timp poate fi de asemenea un interval finit oarecare: 
L“(10, 27]). 


Capitolul 2 - Semnale și Sisteme. Semnale: Clase de Semnale 


Putem “măsura” un semnal ? 
Norme de semnale. Fie u € LI(R), x € (Z). Atunci 


+-oo Too 
hula = J (rar și llh:= X le) 


= k=— oœ 


sunt norme bine definite pe L (R), respectiv pe l! (Z). Aceste norme măsoară “acțiunea” 
semnalului u, respectiv zx. 


In mod similar, rădăcina pătrată a energiei totale a unui semnal u € L? (R), x € l2(Z) 
definește o normă pe L?(R), respectiv pe l2(Z): 


2 


his | [O ratar)) i ll | 5 ewl 


OO k=- 


In fine, norma sup este definită la fel pentru semnalele mărginite, indiferent dacă sunt 
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continuale sau discrete: 


suplu(t)|, supļz(k)| =: zl. 
teR kEZ 

Se observă că în cazul discret norma sup coincide întotdeauna cu norma {°°. In cazul 
continual, este posibil ca norma L% să fie bine definită |ulla := essuplu(t)| < o, 
adică u să fie un semnal esenţial mărginit, dar care să nu fie neapărat mărginit, deci 
norma sup a lui u să nu existe. 


sint dacă teR-N 


Exemplu: u(t) = | t dacă teN 


Avem [lullo = 1, chiar dacă u nu este mărginit. 
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Alte exemple 


Exemplul 5. Fie semnalele continuale 


-3 O<ti<l 
0 în rest 


=i E | 


Este evident că u este un semnal mărginit, cu sup|u(t)| = 1 (și deci cu ||ul|x = 1), dar 
tER 

care nu aparține nici lui L (R), nici lui L? (R). 

v nu este (esențial) mărginit (lim v(t) = +00) și nici un semnal de energie finită, insă 


v € LI(R) şi lloll = 2. 
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Alte masuri 


Media valorii absolute! a unui semnal (continual) este dată de 


In electronică, valoarea unui semnal de curent alternativ este exprimată de rădăcina 
medie pătratică (root mean square=RMS) 


N= 


RM S(u) := [p(u)]?. 


l AA=absolute average 
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De exemplu, pentru semnalul treaptă unitară, AA(1(t)) = 1/2 și RMS(1(t)) = 1/2. 


l l 1 l l l 
Pentru un semnal sinusoidal, RMS este — x amplitudinea semnalului. 


V2 
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Structura spaţiilor de semnale 


e Structura algebrică: spaţii vectoriale. 
e Structura topologică: spații Banach (spaţii vectoriale normate și complete). 


e Structura geometrică: spaţii Hilbert (spații vectoriale normate, complete, cu norma 
definită de un produs scalar). 


Suntem interesaţi de 


- proprietăți calitative: semnal mărginit/nemărginit, convergent către 0 când t — oo, 
periodic, etc. 
- proprietăți cantitative: ||u||, ||u]|2, cât de repede converge la 0, etc. 
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1.3 Operații cu semnale. 


Structură de spațiu vectorial: adunare (f + g), înmulţire cu scalari (a: f, a € R sau C). 
Structură de algebră (Banach): se introduce o operaţie de înmulțire. 


Aceasta poate fi produsul uzual (f x g) sau un alt tip de “inmultire”: de exemplu, 
convoluția. 


Definiţia 6. 


1. Fie u,v € SR. Presupunem că pentru t € R, funcția 7 — u(t—r)u(7) este 
integrabilă pe R. Atunci funcţia 


+00 
w(t) = J u(t — 7r)u(r)dr =: (uxv)(t) (1) 


+o 
D= J u(Ou(t —0)40 = (v * u)(t) 


— OO 


este bine definită în t € R și se numește produsul de convoluție sau CONVOLUTIA 
semnalelor continuale u și v. 
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2. Fie x,y € Si. Presupunem că pentru n € Z, funcția k — x(n — k)y(k) este 
sumabilă pe Z. Atunci funcţia 


z(n) = dU any) =: (£ *y)(n) (2) 


este bine definită în n € Z și se numește produsul de convoluție sau CONVOLUTIA 
semnalelor discrete x și y. 
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Proprietăţi 


1. Dacă u,v € LI(R), atunci u xv este bine definită a.p.t. în R ṣi u*v € LI(R). In 
plus, ||u x ul < llullı llvllı. Se poate arăta că (L! (R), +,-, +) este o algebră Banach. 


2. Dacă he LI(R), u € L?(R), atunci h xu este bine definită a.p.t. în R și 
h xu € L?(R). De asemenea, ||h + ulla < |n lullo. 


Q: Există element neutru la convoluție, adică un semnal d astfel încât f xd = dxf =f, 
pentru orice f ? 


Răspunsul la această întrebare diferită substanțial între convoluția semnalelor continuale 
și a celor discrete. 
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Translaţia în timp 


Operaţie care joacă un rol important în definirea proprietăţii de invarianță în timp 


Definiţia 7. Fier e R (l € Z). Se numește operator de translație (sau shift), operatorul 
07: SR — Sr (o! : S? — SĂ), definit de 


(07u)(t) =u(t-—7), teR  ((olx)(n)=z(n-0), neZ) (3) 
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Figura 13: Translaţie în timp (7 = 1) 
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1.4 Semnale cu Impulsuri (Semnale Singulare) 


e Motivaţie: există situații în care anumite semnale (funcţii) acţionează pe intervale 
foarte scurte de timp, unde pot lua valori extrem de mari. 


e Consecinţă: este imposibil să se măsoare valorile instantanee ale unui astfel de 
semnal (există o limită fizică a măsurării unui interval de timp!); se poate insă 
observa /măsura efectul acțiunii acestui semnal. 


Exemplul 8. Lovirea unei mingi (de tenis, de fotbal). Presupunem ca o forță F 
acționează asupra mingii (de masă m) în intervalul tọ = 2.999sec și tı = 2.001sec. 
Presupunem că la momentul to mingea se afla în repaus. Efectul acestei acţiuni este dat 


de 


3.001 3.001 
| F(r)dr = m v(3.001) — m u(2.999),  v(3.001) = v(2.999) + Pdr 
2.999 M J2.999 


Cu alte cuvinte, putem observa efectul forței F pe intervalul considerat, măsurând viteza 


mingii după încetarea acțiunii forței F. 
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Exemplul 9. Incărcarea rapidă (instantanee) a unui condensator. Se consideră circuitul 
din figura de mai jos (figura 14). 


N t=0 


i(t) 
1V E v(t) | C=1F 


Figura 14: Incărcarea unui condensator 


Presupunem că v(0) = 0 și q(0) = 0. La închiderea circuitului tensiunea la bornele 
condensatorului C, v(t), crește (aproape instantaneu) la valoarea V, iar curentul în circuit 
va atinge o valoare foarte mare, după care va fi practic nul (vezi graficele din figura 16, 
când R — 0). De asemenea, sarcina va fi transferată la bornele condensatorului aproape 
instantaneu, 


tot = i iri, =, daca C=1F; V = IV. 
0 
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Se constată că deși i este 0 a.p.t (i(t) = 0 pentru t Æ 0, i(t) ~ œ pentru t = 0), avem 


| ioa 70! 


Am neglijat cu totul rezistența existentă în circuit. 


Concluzii. Astfel de semnale (F, dar mai ales i) au un comportament impulsiv. Ele 
se mai numesc și semnale impulsive sau singulare. Studiul semnalelor singulare este 
încadrat de Teoria Distribuţiilor (L. Schwartz, 1950). 


Ne vom limita la studiul (din punctul de vedere al ingineriei electrice, și nu al teoriei 
distribuțiilor) unui singur semnal singular, impulsul Dirac. 
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Impuls Dirac 
Definiţia 10. Se numește impuls Dirac, notat ô(t), (un “obiect” care este) o idealizare 
a unui semnal având proprietățile: 


a. este foarte mare intr-o vecinătate a lui t = 0: 6(t) este nedefinit în O; poate fi chiar 
infinit. 


b. este foarte mic în afara acestei vecinătăți: 6(t) = 0 pentru t 70. 


C. [o ad = d; 


— OO 


Capitolul 2 - Semnale și Sisteme. Semnale: Semnale Singulare 


Circuitul RC 


Notăm cu R rezistenţa în circuitul din figura 14. 


Figura 15: Circuit RC 
DI i Ri = V — v rezultă 


dv _ 1 1 
dt RC RC 
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dacă v(0) = 0, soluţia acestei ecuații se poate scrie 


Eon 


v(t) = V(1 — e-e) T 


-2 
l—e R. 


Figura 16: Tensiunea și curentul condensatorului 
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Se constată că lim v(t) =); imi = 


R—0 


0 170 
œo t=0 ` 


Este 6(t) = lim i(t) ? NU, deoarece 


—0 


Pe de altă parte, 
-+- OO 
lim (/ măr) = lim (CV) =, 


R—O 23 


adică ip(t) este o “aproximaţie” a impulsului Dirac 6(t) pentru valori “mici” ale lui R 
(în condiţiile în care C = 1, V = 1). 


Capitolul 2 - Semnale și Sisteme. Semnale: Semnale Singulare 


Aproximaţii ale lui ô 


In figura 17, avem p-(t) = 1(t) — 1(t— £), de unde 


Se poate arăta riguros că, în sens distribuțional, impulsul Dirac ô(t) este intr-adevăr 
derivata treptei unitare 1(). 


Nu contează forma și valorile pe care le ia o aproximaţie oarecare a lui ô, ci efectul 
acțiunii acesteia, adică faptul că | = 1. 
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P (4) 


1/75 


52 ep t 


Figura 17: Aproximații ale lui ô 
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Mai precis, distribuția (impulsul) Dirac este o funcțională (pe spaţiul funcţiilor test) 
definită prin 
FOO 


wm= FORMAL J OE de: (4) 


= 
Se poate defini în mod similar distribuția asociată semnalului treaptă unitară, 


10) | led | rar 


— OO 


Să mai notăm că, în conformitate cu (4), impulsul Dirac este element neutru la convoluţie 
pentru semnalele continuale (se ia f(T) := h(t — 7), t fixat). 


Spre deosebire de ô (semnal singular), elementul neutru al operației de convoluţie discretă 
este impulsul discret, care este un semnal regulat. 
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2. Transformari Integrale 


Transformarea Fourier + Fourier discretă 


Transformarea Laplace + Laplace discretă sau transformarea Z 


Transformarea Fourier. 
Definiţia 11. Fie f € L!(R), f contină. Atunci funcţia 


F:jR—C, F(ju):= i Joe mų a (5) 


este bine definită, continuă și mărginită pe R, și se numește transformata Fourier a lui 
f în punctul jw. 

Aplicația f — F, F = F(f), se numește transformarea Fourier; F este un operator 
liniar, F : LI(R) — CO, unde C? este mulțimea funcțiilor continue și mărginite pe R, 
având limitele la 4+oo egale cu 0. 
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Interpretare fizică 


Notă. Am notat domeniul de definiţie a lui F cu jR (în loc de R), respectiv argumentul 
funcţiei F cu jw (în loc de w) pentru a sublinia diferenţa între axa reală a momentelor 
de timp și axa reală a frecvenţelor. 


Spectrul în frecvenţă al semnalului f(t) 


jarg| F'(ju)] 
F(ju) = |FUu)| e tz (6) 


amplitudinea 


Amplitudinea (în frecvență) a lui f(t) se măsoară în decibeli (dB) iar faza (în frecvență) 
a lui f(t) se măsoară în radiani. 


Formula (5) poate fi interpretată ca o combinație liniară de oscilații armonice e'“t de 
amplitudine variabilă |F (jw)|. 


F este o rezoluţie de frecvență a lui f, evidențiind amplitudinile (|F'(jw)]) oscilaţiilor 
armonice (e'“!) din care este compusă f. 
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Proprietăţi 


+00 
0. Formula de inversare. Dacă F e L!(jR), adică J F(ju)| dw < œœ, atunci 
transformata Fourier inversă a lui F, f = F-1(F), este bihe definită și dată de 
BE | 
F:R-—C, ft) ==] F(jw) e%** dw. (7) 
T — OO 


1. Liniaritate. Faf + pg =aF(f)+6F(g). 

2. Scalarea axei de timp (asemănare).  f(at) A 1 F), a > 0. 
3. Translaţie în timp.  f(t—7) Led F(jw), TER. 

4. Translație în frecvență. ef^ f(t) 1z F(j(w —A)), AER. 


5. Convoluţie în domeniul timp. (Produs în domeniul frecvență). 


fxg > F(f) Flg) =FG. 
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(F=) = FFG) => 


+00 
i . jwt 
=) F(ju)G(ju) e“ dw. 


6. Produs în domeniul timp. (Convoluţie în domeniul frecvență). 
J= OO f , 
fg > =z; FU — Ga) da. 


7. Egalitatea lui Parseval. Fie u € LI(R)N L2(R) (care este o mulțime densă în L2(R)) 
și U = F (u). Atunci 


1 l ] -+- oo | r 3 
ul = -lU Gole = Za | i Ie) d) 


iar transformarea u > Jaf (u) este un izomorfism de spaţii Hilbert, de la L2(R) la 
L?(jR). 
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Convoluţie. Filtrare. 


Ecuațiile y(t) = (h x u)(t) Y (ju) = H(ju) U(ju) descriu un filtru de frecvenţă. 
De exemplu, sistemul auditiv se comportă ca un filtru. 


Convoluţia “distruge” toate frecvențele (oscilaţiile armonice) care intră în componența 
lui u, dar care nu apar în h. Sunt frecvențele la care H(jw) este O sau foarte apropiat 
de 0. 


Filtrare & predicţie: N. Wiener 1930. 


Fie h € LI(R), cu H(ju) 70, Vw e R și fie de asemenea un semnal y € R. Atunci, 
+00 


pentru orice e > 0, există u cu proprietatea că J (h x u)(t)— y(t)ldt < e. 


— OO 
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Transformarea Fourier discretă 


Considerăm mulțimea de momente de timp finită 7 = {0,1,2,..., N— 1}, N >2. 
Orice șir (x, )nez se numește semnal finit cu N eșantioane. 


Definiţia 12. Se numește transformare Fourier discretă a semnalului (£n)nez un alt 
semnal finit (Xk)kez definit de 


N-—1 N—1 
27 __j27 
As X Tne INEN — X pa, aqi=e ÎN, (8) 
n=0 n=0 


XX se numește eșantionul spectrului lui x pe frecvența k. 
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Dacă notăm z” = |zo zi ©- 2-a] și XI := [Xo Xa ©- X-a] atunci 
relația (8) se poate rescrie în formă matriceală 


1 1 1 

1 a g~ 
X=Wzr, unde W := |1 a? aa N= 

| i a aN -1)x(N-1) 


Deoarece W W = W W = N Iy, rezultă imediat formula de inversare a transformatei 


Fourier discrete: i 
= —W X. 
“=N 


Cu W s-a notat conjugata matricii W. 


Transformata Fourier discretă este utilizată la calculul transformatei Fourier. 
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Presupunem că x(t) este nulă în afara unui interval J. Se aleg (in 1) N eșantioane 
(N = 2P) și se obține (£n)nez. Se calculează eșantioanele spectrului X (jw) în punctele 
0, 27/N, 4r/N, ..., 27 (N — 1)/N, 


2 27 
XGL Rk) = J r(t) e INF dt, 


N I 


Integrala se aproximează cu ajutorul unei sume în care apar cele N eșantioane ale lui x. 


Transformata Fourier Rapidă: Procedură de calcul eficient. 


Detalii complete la cursul de Prelucrare Numerică a Semnalelor. 
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Transformarea Laplace 


Consierăm transformarea Laplace unilaterală (la dreapta). 


Definiţia 13. Fie f € Sg. Se numește transformata Laplace unilaterală la dreapta a lui 
f în punctul s 


F(3) := J aed (9) 


F este bine definită în s (integrala improprie converge) dacă s € Sf, 
unde S% = {s € C : (9) este absolut convergentă }. 


i i L ; F 
Aplicația f =} F se numește transformarea Laplace (unilaterală la dreapta). 


Notă: Se pot defini în mod similar transformatele Laplace unilaterală la stânga (£L_) și 
bilaterală (£): orizontul de integrare se ia în (9) (—o0,0), respectiv (—o0, 00). 


Notaţie: Vom nota de aici înainte pe L4 cu £; discutăm numai despre transformarea 
Laplace unilaterală la dreapta. 
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Funcţii original (Laplace) 


Introducem clasa funcțiilor original. Spunem că f : R — C este o funcție original 
(Laplace), f € O, dacă f are următoarele proprietăţi: 


(i) f(t) = 0, pentru t < 0. 
(ii) f este continuă pe porțiuni în [0, o0). 


(iii) există M > 0, so > 0 astfel încât |f (t)| < M e**, pentru orice t > 0. Numărul 
real sg se numește indice de creștere (exponențială). 


Teorema 14. Fie f€ O o funcție fixată de indice so. 
Atunci S; = {s € C : Res > so} şi F(s) este olomorfă în S$. 
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Notaţie. Terminologie. Transformata Laplace a funcţiei f în punctul s: F(s) = 
LI f(t))(s). Transformarea Laplace: F = L(f). 
Funcţia F se numește funcţia imagine (Laplace) a funcţiei (original) f. 


Exemplu. Funcţia treaptă unitarăl(t) este o funcţie original având so = € > 0. 
Atunci Ss? = {s €C : Res >Q} şi 


L{10)}(6) = J "etat = L 


S 
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Proprietăţi 
0. Transformata Laplace inversă. Fie f € O cu indicele de creștere sp și F = £(f). 
Atunci | 
1 +70 
T(î) = = | F(s)e*! ds, Vo > so, şit>0. (10) 
21] 0—]00 


1. Liniaritate. £Liaf+bgi(=)aFr+pbG. 


2. Asemănare (scalarea axei de timp). £Lif(at)!(s) = Er ), 


S 
Q 


3. Translaţie (întârziere) în timp. L{f(t—rT)}(s) =e 75 F(s), T >0. 
4. Translaţie în frecvență. £Lle-“ f(t) (s) = F(s +a). 

5. Derivarea imaginii. L{t” f(t))(s) = (—1)" F™ (s). 

6. Derivarea funcţiei original. Dacă f, f, „fim e O, atunci 


LIf(0)(3) = s"P(3) — 8%" O —... — F0%9(0), 
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unde g(04) = lim g(t) este limita la dreapta în 0 a funcţiei g. 
In particular, L{ f (t)}(s) = sF (s) — f(0.). 


7. Teorema valorii finale. Dacă f,f € O și dacă există lim F(t) not Flo) w, 
atunci lim sF(s) = f(o0). 


8. Teorema valorii iniţiale. Dacă f, f € O și dacă există lim sF(s), atunci 
lim sF (s) = f(0.). 


9. Convoluţia. Dacă h,u € O și h xu este bine definită, atunci y = h x u este funcție 
original și 


Liylt) ils) = H(s)U (s) (11) 
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1. e(t1(5)(s) £ a Cea LANs) È a 
2. L{e* 1(5))(s) £ = 
3, 
PT —jwt ıl iwt l —jwt 
tcosut)(s) = LE + eI} (8) E Eeeh) + ce” its) 


S O a 
2s—ju  2s+ijw a 4u 


I> 


In mod similar se obține £isinuwt!(s) = pier 
s*a 
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Transformata Z (Laplace discretă) 


Definiţia 15. Fie x : Z — R(C) un semnal discret. Se numeste transformata Z 
unilaterală la dreapta a lui z în punctul z, functia 


e e lat (12) 


k==00 
definită pe Si - domeniul de convergenta al seriei (12). 
Vom nota in mod obisnuit X(2) = Z{æ(k)}(z), sau, mai simplu, X = Z(x). 
Aplicatia x -— X se numeste transformarea Z. 
Domeniul de convergență S1 este in mod uzual o coroană circulară. 


Detalii complete: TSII - capitolul Sisteme discrete. 


Capitolul 2 - Semnale și Sisteme. Transformari Integrale 


3. Sisteme. Proprietati Fundamentale. 


Definiţia 16. Inţelegem prin sistem în accepțiunea intrare-ieșire o aplicație T : U —> yV, 
y = T(u), unde U, Y sunt spaţiul semnalelor de intrare, respectiv spațiul semnalelor de 
ieșire. 


Sistemul se numește (cu timp) continuu dacă U și Y sunt spaţii de semnale continuale, 
(cu timp) discret dacă U și Y sunt spaţii de semnale discrete, respectiv hibrid dacă U și 
Y sunt unul spaţiu de semnale continuale iar celălalt spaţiu de semnale discrete. 


In mod uzual U, Y € (Sp, Sc, S8, Să). 
Caz remarcabil: U = L2(R) (l2(2)), Y = L2(R) (2(2)). 
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Sisteme liniare 


Definiţia 17. Fie sistemul T : U — YV, y = Tu, unde U și Y sunt spaţii liniare. Dacă T 
satisface principiul superpozitiei, 


T (œu + azua) = aiT (u1) + œT (u2), V u1, u2 E U, V ali, Q2 € R(C), 


atunci T este un sistem liniar. 


Exercițiu: arătaţi că circuitul RC definește un sistem liniar. 
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Invarianţă în timp 
Definiţia 18. Un sistem T : U —Y, y = T(u), este invariant în timp dacă 


o'T ="To',vreR, (13) 


unde o” este operatorul de translație (shift) introdus în (3). 


Această proprietate este ilustrată în figura 18 care poate fi interpretată astfel: dacă 
la intrarea u (arbitrar alesă dar fixată) răspunsul sistemului este y = Tu, atunci se 
cunoaște răspunsul sistemului la orice intrare de tipul u(t—7), 7 € R: acesta este 


y(t —7). 


Intr-adevăr, din (13) rezultă că o“Tu = Tou, sau echivalent, oy = y(t — T) = 
T(u(t—7)). 
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Figura 18: Invarianţă în timp 
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Exemple 


1. Fie sistemul y = T(u), y(t) = u?(t). 

Este sistemul liniar ? Nu. De exemplu, T(3 u) = 32u? Z 3T(u) = 3u?. 

Este sistemul invariant în timp ? Da. Fie u un semnal de intrare oarecare și fie 
y(t) = u2(t) semnalul de ieșire corespunzător. Fie r € R; notăm u(t):=u(t-—7), 
y (t) := y(t — T). Atunci T(u) = u? (t — T) = y(t — T) = yr, aşadar sistemul este 
invariant în timp. 


2. Fie sistemul y = T (u), y(t) = t? u(t). 

Este sistemul liniar ? Da. [T (au+ 8v) (t) = t (au(t)+8v(t)) = al (uN (AÐ + 8T (v) (t). 
Este sistemul invariant în timp ? Nu. Considerăm u(t) = 1 (t). Atunci y(t) = t°1(t) și 
y(t) = (t = T) L(t — T); insă y(t) = [T (u-)](t) = P10- 7) # y(t). 
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Cauzalitate 
Definiţia 19. Un sistem T : U — V, y = T(u), este cauzal dacă: 


1. ieșirea la momentul 7, (7), nu depinde de intrarea sistemului u(0) evaluată la 
momente de timp ulterioare lui 7, 6 > 7. 


sau, echivalent 
2. u(t) = uo(t) pentru orice t < 7, implică yi(t) = y2(t) pentru orice t < 7. 
sau, echivalent 
Sa dpi pu E pal. O dal lee = 0, pentru oer ER. 
Operatorii de la punctul 3. sunt familiile de proiectori P, +, T € R, 


ut) tT, SA. 
(Pr) = e Palos ed 
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y(t) 


t (Ru XÐ 


Z > 


Figura 19: Cauzalitate 
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Observaţia 20. Dacă în definițiile 18 și 19 înlocuim pet cu k, per cun și respectiv pe 
R cu Z, obţinem definițiile invarianţei în timp și cauzalităţii pentru sistemele discrete. 


Exercițiu. Este sistemul discret y(n) = u(n — 1) + u(n — 1) cauzal ? Dar liniar și 
invariant în timp ? 
Scrieţi definițiile 18 și 19 pentru sisteme discrete. 


Se poate “descrie” un sistem care este liniar, invariant în timp și cauzal ? 
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4. Sisteme de Convolutie 


In anumite situații (care depind de proprietăți ale spaţiilor de semnale U și Y) se poate 


arăta că un sistem liniar și invariant în timp este un sistem de convoluție. 


Un sistem y = Tu este un sistem de convoluţie dacă există un semnal Ph astfel încât 
y = hxu. Funcţia h se numește funcţia pondere a sistemului de convoluţie. 


Sistem de convoluţie cu timp continuu: 


+00 Too 
TTE i îi jul E ] (Out — 0)48. 


OO — OO 


Sistem de convoluţie cu timp discret: 


+00 +00 
y(n) = 5 h(n-ku(k)= 5 hun- 1). 
k=— o0 l=—0 
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Sisteme de Convolutie 


Propoziția 21. Un sistem de convoluţie este liniar şi invariant în timp. 


Demonstraţie.  Liniaritate. Dacă yi(t) = (h x u1) (t) și y2(t) = (h x uo)(t), atunci 


+00 
j J E E 


— OO 


[o h(t — Thu (7)d7 + [o h(t — T)uo(7)d7 


— OO — OO 


(h x u1)(t) + (h * u2)(t) = pat) + y2(t). 


Invarianță în timp. Fie u(k) o intrare oarecare a sistemului de convoluție discret 
y(n) = (hxu)(n) și fie ğ(n) ieșirea sistemului la intrarea ŭ(k) = u(k— 1), l € Z arbitrar, 
fixat. 
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Avem 


-+-oo T 
In) = X hn-kūlk)= X h(n- kūlk-—l) 
k=—00 k=—00 
k —l= mata 
= XO h(n-l-m)ŭ(m) = y(n- l) 


Aşa cum am menționat mai înainte, reciproca acestei propoziții nu este, în general, 
adevărată. 


Caind este un sistem de convoluție și cauzal ? 


Propoziția 22. Un sistem de convoluție y = h x u este cauzal dacă şi numai dacă 
h(t) = 0 pentru oricet <0 (h(n) =0 pentru orice n < 0). 
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Răspuns în timp al sistemelor de convoluţie 


A. Răspunsul la impuls. 


1. Cazul discret. Fie y(n) = (h x u)(n) un sistem de convoluție discret. Dacă 
u(k) = 6(k) (impulsul discret a fost introdus în secțiunea 2.1.1) atunci răspunsul 
sistemului la impuls este 


+00 
yi(n)= X h(n- k)ô(k) = h(n). (14) 


k=— oo 


2. Cazul continuu. Fie y(t) = (h x u)(t) un sistem de convoluție cu timp continuu. 
Dacă u(7) = 6(7) (impulsul Dirac a fost introdus în secțiunea 2.1.4) atunci răspunsul 
sistemului la impuls este 


+00 
Ei J h(t — r)â(r) dr = h(t). (15) 
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Semnificaţia funcţiei pondere 


Notă. Justificarea celei de-a doua egalităţi din (15) se face în cadrul teoriei distribuţiilor. 
Se ia f(T) = h(t — 7) în egalitatea (4) și se obține (15) (vezi și secțiunea 2.1.4). 


Funcţia pondere a unui sistem de convoluţie este răspunsul la impuls al 
sistemului respectiv. 


Capitolul 2 - Semnale și Sisteme. Sisteme de Convolutie 


B. Răspunsul la treaptă unitară. 


1. Cazul discret. Fie y(n) = (h x u)(n) un sistem de convoluţie discret, unde 
u(k) = 1(k). Atunci răspunsul sistemului este 


p(n) = X h(n-k)i(k)= X R) (16) 


k=- [== 


2. Cazul continuu. Fie y(t) = (h x u)(t) un sistem de convoluție cu timp continuu, 
unde u(7) = 1(7). Atunci răspunsul sistemului este 


Ap 00 i 
yı(t) = J h(t — T)1 (T) dr = J h(0) dð. (17) 


— oo — 00 


— dyı(t) 


Se observă că h(t) = yı (t) d 


treaptă unitară. 


: funcţia pondere este derivata a răspunsului la 
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5. Raspuns in Frecventa. Functie de Transfer”. 


C. Răspunsul la intrare de tip armonic. Fie u(t) = e*t, unde w = 27 f este pulsația iar 
f este frecvența semnalului armonic. 


FOO +00 
Yalt) J h(0)eie(t—6) dô = ] h(0)e- io? do eiet 


— OO — OO 


h(jw) u(t). (18) 


Funcția ȘI 
(ju) = J h(0)e 7i“? 40 


— OO 


este răspunsul în frecvență al sistemului de convoluție y = A x u. 


? Incepând cu secțiunea 2.5 ne vom concentra atenția în exclusivitate asupra sistemelor cu timp continuu. Cazul discret 
va fi reluat în capitolul dedicat sistemelor discrete. 
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Proprietăţi 
Observaţia 23. 


1. h(jw) = H (jw), unde H (jw) este transformata Fourier a funcției pondere h. Aceasta 
este bine definită dacă h € L! (R). 


2. yn(t) este un semnal armonic de acceaşi frecvență w cu intrarea u(t), dar de 
amplitudine și fază modificate. 


3. Mai precis, fie u(t) = A coswt şi h(jw) = |h(jw)| efr IRG), 
Atunci se poate arăta că (exercițiu) 


y(t) = Afh(jw)| cos(wt + argfh(jw)]). (19) 


4. Dacă h(t) ia valori reale, atunci se poate demonstra că (exercițiu) 


h(jo) = h(-ju), 
de unde |h(juw)| = |h(—jo)|, arglh(ojo)] = —arglh(ju)]. 
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D'ivizorul de tensiune: reluat 


Figura 20: Divizorul de tensiune în curent alternativ 


Exemplul 24. Ecuațiile (cu u = vi, y = vo) 


di R. 1 | 
di = S + AL i(0) = 0 
y = du 


t 
1 
Rezultă i(t) = i T i zu(T)dr, de unde 
O 
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Egalitatea de mai sus arată că circuitul definește un sistem de convoluție (cauzal). 


| R | 
Funcţia pondere este h(t) = rată pentru t > 0 și A(t) = 0 pentru t < 0. 
Funcția de răspuns în frecvență este 


p PR R R 
(iu = aa e N AA 
iu) = | gat BLER 


ceea ce arată că circuitul LR are un comportament de tip filtru trece-jos. Intr-adevăr, 
relația (19) arată că amplitudinea ieșirii vọ(t) la o intrare de tip armonic de pulsație 
“mare” (de exemplu, v;(t) = A cosuwt cu w > 10*R/L) este practic 0. Filtrul “distruge” 
amplitudinile oscilaţiilor armonice de frecvență mare, pentru care h(ju) = 0. 


R 


— ———V;(s) caracterizează comportamentul I/O în 
sL+ R il ) P / 


Pe de altăparte, relația Vo(s) 
domeniul operaţional. 


H (s) : 


R , i KEET. 
:= ——— este așa-numita funcție de transfer a circuitului. 
sL+ R 
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Funcție de transfer 


Fie sistemul de convoluţie y = h x u. Fie s € C și u(t) = e. Răspunsul sistemului este 


+00 +00 
y(t) J h(O)estt-0) dð = ] h(0)e- 5 d0 es 


— OO 


| 
= 
(Va) 
pl, 
= 
P= 
œ 
s 


(20) 


Funcția 


+00 

H(s) = J h(te-** dt (21) 
— OO 

este funcția de transfer a sistemului de convoluție y = h * u. 

Aceasta este bine definită în punctele s € C pentu care integrala din (21) converge 


absolut, adică în S}. 


Dacă sistemul de convoluție este cauzal (A(t) = 0 pentru t < 0), atunci 


H(s) = | h(te-** dt (22) 
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este transformata Laplace (unilaterală, la dreapta) a funcţiei pondere h(t); 
domeniul de definiţie: S. 


Dacă S} c (S, ASY) și U(s) = L{u(t)}s), Y (s) = L{y(t)}s), atunci 


Domeniul frecvență 


Domeniul timp 
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Observaţii 


1. In mod uzual, vom considera sisteme de convoluţie cauzale ale căror funcţii pon- 
dere sunt funcţii original (Laplace). De pildă, ecuaţiile diferențiale liniare cu coeficienţi 
constanţi definesc astfel de sisteme (vezi exemplul ). Cum semnalele (de intrare) standard 
sunt de asemenea funcţii original, u € O, rezultă că y = hxu € O și Y (s) = H(s)U (s). 


2. De altfel, egalităţile (23) sunt valabile într-un context mai larg, impus de situaţia în 
care sistemul este definit de o ecuaţie diferențială. Vom considera în cele ce urmează 
sisteme de convoluţie având funcţii de transfer raţionale, definite de ecuaţii diferenţiale 
liniare având coeficienţi constanți. 


3. Formulele (23) sugerează o metodă de calcul (analitic) a răspunsului unui sistem 
pentru o intrare dată. Cum H(s) și U(s) se cunosc, Y (s) rezultă în mod banal și y(t) 
se obține ca transformată Laplace inversă a lui Y(s). 
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0. Sisteme de Convolutie cu Functii de Transfer Rationale 


e Ecuațiile fizico-matematice ale unui sistem: exprimarea matematică a relațiilor dintre 
variabilele care intervin în sistemul fizic. 


e Relaţiile dintre variabile --> ecuaţii diferenţiale: exprimă în mod obișnuit o ecuaţie 
de echilibru, determinată de principiile (legile) fizicii care descriu fenomenele care au loc 
în sistem. 


e Sisteme fizice cu parametri constanţi. Elemente liniare, intervale de liniaritate, 
liniarizare. Sisteme fizice cu parametri concentrați. 


Modele matematice: 


Ecuaţii diferenţiale ordinare, liniare și cu coeficienţi constanți. 
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Exemplu: circuitul RLC serie 


- intervin elemente liniare de circuit: R, L, C (de exemplu, R este un element liniar 
numai în intervalul de funcţionare pentru care a fost prevăzut: legea lui Ohm). 


- coeficienţii sunt constanți (ipoteză de lucru: capacitatea condensatorului variază foarte 
lent cu timpul, fiind practic constantă). 


Relația dintre variabile (scrisă în baza legilor lui Kirchoff și a caracteristicilor curent- 
tensiune la bornele elementelor R, L,C): ecuaţie diferențială ordinară, liniară și cu 
coeficienți constanți. 


d2vo dvo 
“pRO 


O le 
Ss di? dt 


+ Vvo = v(t). 


In accepțiunea I/O: 


v; (tensiunea la bornele circuitului) - intrarea sistemului 
Uo (tensiunea la bornele condensatorului) - ieșirea sistemului. 
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Ecuația (forțată sau comandată) definește (în mod riguros) un sistem. 


In ce manieră ? 
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Ecuaţii diferenţiale (forțate) 


Se consideră ecuaţia diferenţială în y: 


i E a ee ae 
nn | În gpn=i din “O 
daul) (d) du(t) 
= m i= Ll xa a e.. 6] Ên $ 24 
b Jem b 1 dim= T bi T bou(t) a Z 0 ( ) 


Numărul n € N se numește ordinul ecuației diferențiale. 


Rescriere compactă (D este operatorul diferențial): 


A(D) yt) = B(D) ut), (Dx)(t) = 


Aici A(s) = ans” +... + ao, B(s) = bms™ +.. . + bo sunt polinoame cu coeficienți reali 
(sau complecși - foarte rar!). 
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“Reţeta” clasică de rezolvare 


1. Se determină soluția ecuației omogene A(D) y(t) = 0 având ecuația caracteristuică 
A(s) = 0. Fie A, ... , A, rădăcinile acesteia de multiplicități ma, ... , My. 


Soluţia ecuației omogene are forma generală 
rT 
A;t 
Volt) = N milt) €e ; 
i=1 


unde 7;(t) sunt polinoame de grad cel mult m; — 1 în t, ai căror coeficienți se determină 
din condițiile inițiale. 


2. Se determină (pentru u(t) specificat) o soluție particulară yp(t) a ecuației complete 
A(D) y(t) = BID) u(t). 


3. Soluția generală este 
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Rezolvarea ecuaţiei diferenţiale cu ajutorul transformării Laplace 


Se aplică transformarea Laplace în ambii membri ai ecuației (24). Cu notațţiile obișnuite 
Llu(t))(s) = U(s), Liyl(t)i(s) = Y (s), se obţine (utilizând proprietatea 6. a trans- 
formării Laplace) 


an[s Y (5) — s”-ty(04) = y- (04)] +... + cas (8) — y(0+)] + aY (8) 
= ba Ues ula mat O O Oa] +... + ba[sU(s) — ala + boU(s) 


sau, în formă compactă, 


unde AKo) =a (0 H lany O) taan Oa a Pan VOnt 
a1y(0+)), Bi(3) = bmu(04)s®1+[bmu' (0+)+am-18(04)]57?+. [bul Y (04)+ 
... + bru(0.)] sunt polinoame de grad cel mult egal cu n — 1, respectiv m — 1. 
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Rezultă 


B(s) A;(s) B;(s) 
Y = U — 25 
= aj 9 aţi) Alo) = 
Componentă forțată  Componentă liberă Condiții inițiale u 
In conditii initiale nule, y(0,) =y (0) = ... =g a0 =0 a0 = w(0,) = 


„= uY) (04) = 0, ecuaţia (25) devine 


Y(s) = H(s)U(s), unde H(s) := 


este funcția de transfer definită de ecuaţia (24). 


De cele mai multe ori, în mebrul drept al ecuaţiei diferenţiale (24) nu apar derivatele lui 
u (m = 0), situaţie în care (B(s) = bo, B;(s) = 0) 
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Consideraţii finale 


“Morala”: Ecuația diferențială (24) definește un sistem de convoluţie, y = h x u, având 
B(s) 


„a cărui funcție pondere este dată de 
A(s) 


funcție de transfer raţională, H (s) = 


h(t) = L~ (H (s))(t). 


Acest sistem u — y realizează tranziția între semnalul de intrare (comanda) u și 
componenta forțată a răspunsului sistemului (soluției ecuației diferențiale). 
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CAPITOLUL 3: SISTEME CONTINUE SISO 


PREAMBUL: Funcţii raţionale O funcţie raţională 


ms +... +bis+bo _ B(s) 


H (s) = — 
(5) ans” +... +ais+ao A(s) 


se numeşte: 


- (strict) proprie, dacă (m <n) m<n. 
- biproprie, dacă m =n. 
- improprie, dacă m > n. 


Presupunem că H(s) este o rațională ireductibilă. 
Rădăcinile polinomului B(s) sunt zerourile finite ale lui H(s). 
Rădăcinile polinomului A(s) sunt polii finiți ai lui H(s). 


Zerourile/polii la infinit ale raționalei FI (s) sunt zerourile/polii în O ale raționalei H (2) 
S 
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Se consideră că o astfel de rațională are un numar egal de poli și zerouri (numărând 
multiplicitățile și zerourile/polii de la infinit). 
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| o. | 25 — 1 iad FE 
De exemplu, o rațională strict proprie, H(s) = PETE: are 3 poli finiţi, un zero finit și 
S 
două zerouri infinite. 
| e ji (2 — s)s? p A 
Se vede imediat că H | — | = —— are două zerouri în 0. 
S 1 + s 


O raţională strict proprie are n poli finiţi, m zerouri finite și n — m zerouri la infinit. 


Raţionalele biproprii au toţi polii și toate zerourile finite, în timp ce raționalele improprii 
au n poli finiţi, m — n poli la infinit și m zerouri finite. 


Considerăm în cele ce urmează sisteme de convoluţie cauzale care admit funcţii de 
transfer raţionale, proprii, având coeficienţi reali. 


Acestea sunt sisteme cu o intrare și o ieșire: sisteme SISO (Single Input/Single Output). 


Acest capitol este dedicat ANALIZEI sistemelor continue cu o intrare și o ieșire. 


Capitolul 3 - Sisteme continue SISO Preambul 110 


1. Raspunsul sistemelor SISO 


SISTEM = H(s) - rațională proprie cu coeficienţi reali 
y = hxu Y(s) = H(s)U (s) 
Se dau h și u - H(s) și U(s). Se cere y(t) = L(Y (s))(t). 
Y(s) este o funcţie rațională strict proprie cu coeficienţi reali. Cum se calculează 


transformata Laplace inversă a unei astfel de funcţii ? 


Problema. Se dă o raţională strict proprie cu coeficienţi reali 


F(s) = Bms” +... Past Po _ r(s) 
AnS” +... +s +ao pls) 


Se cere 
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Descompunere în fracţii simple: Teorema dezvoltării. 


Fie p; polii lui F(s) (rădăcinile lui p(s)). Distingem mai multe cazuri: 


I. Poli reali simpli. 


ÎN PNI + CAINII. NE 8 An 
die Ga 5 — S aaa (8 Dr) 5S—P1 i S — p2 id S — Pn (26) 
unde 
A; = (5 = pi) P(5)lozpi TO (27) 
| | i p'(pi) 
Rezultă 


TO= CT — Aa T + Ap eP2 +... + A ePrt. 
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||. Poli simpli, reali sau complecși. Polii complecși apar în perechi complex conjugate. 
Fie o astfel de pereche p1 2 = o £ jw, po = Pi. 


Deoarece toți polii sunt simpli, descompunerea (26) este în continuare valabilă, 


A A A An 
a o ela E mut ea 
SI Sa 93 


o= Pi 
SE r(p2) __ r(Di) = 
Conform relaţiei (27), A> = = = Au. 
Gil i, p'(p2)  p'(Pı) i 
Rezultă 
[ol ( A] 4+ A> (t) — A ePit E Pai — e7*( Ajet“ A je 17) 
S— Pı 5 —p2 


= 2e%™%|Aı| cos(wt + 1), fa = arg[ Aa]. 
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Raspunsul Sistemelor SISO 


Dacă F(s) are 2m poli complecși simpli (adică m perechi complex conjugate, 
Poi—12i = Ci E jwi, i = 1 : m) și n — 2m poli reali simpli, atunci 


F(s)=y, a + 


im T 5 
{=l Pi Pi i=2m+1 


Aşadar 
= ` 2e%i| A;| cos(w;t + arg[44]) + > A, eP”. 


i=l i=2m+1 


II. Poli multipli. 
Fie p1,p2,...px polii lui F(s), reali sau complecși, avân d multiplicitățile m, M2, .. . Mk, 
Mi +t mia +F tr PL Nn: 


Descompunerea lui F se face astfel 


A O II Au 
F(s) = en A n 


Gals =P s= p2) 
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unde 


Se obţine 
k Mli 4—1 k 
dA d B Ail 13 ePit — N milt) ePit. cu grad|r;] < mi. 
= l {=i 


Notă. Am coniderat subînțeles faptul că f(t) = 0 pentru t < 0. Altfel spus, 
A W 
L s Aa £ (o 
5 0), ca 
punerea excesiv, am suprimat termenul 1(t) în toate expresiile de mai sus. 


) = sinwt 1(t), etc. Pentru a nu încărca ex- 
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Analiză calitativă 
Fiecare pol al lui F(s) contribuie cu unul sau mai mulți termeni la f(t). 


- un pol în s =a. Termen de tipul e“. 
v t— 00 . ú t— 00 
Dacă as 0 e” S oaia dacă az 0, e S 0. 


- o pereche de poli complex conjugați în s = o + jw. Termen de tipul e” cos(wt + $). 
Oscilație amortizată (o < 0), cu amplitudine crescătoare (o > 0) sau cu amplitudine 
constantă (o = 0). 


- un pol de ordin m în origine s” = 0. Termen de tipul tm. 


- un pol de ordin m în s = a. Termen de tipul tm Lea. 
v z t— 00 i X Z t—> oo 
Dada- 0T e "= aoar daca < 0r ee 0. 


- o pereche de poli complex conjugați în de ordin m în s = o + jw . Termen de tipul 
e7tt™— lejot, Oscilaţie cu amplitudine descrescătoare (o < 0) sau crescătoare (o > 0). 
Pentru o = 0 nu apar oscilații de amplitudine constantă, ci de amplitudine crescătoare. 
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Ce se întâmplă când F(s) nu este strict proprie ? 
Trebuie extrasă mai întâ partea polinomială, prin împărţirea lui r(s) la p(s): r = 


Avem 


F(s) = q(3) + F(s) = q(s) + a2 cu F(s) strict proprie. 


p(s) 


R? 


Deoarece F(s) este (în situațiile care ne interesează) o rațională proprie, rezultă 
gls) = Omana = E (oo) = const. 


In concluzie, în cazul m = n, avem 


L (F(8))(t) = F(00)8 (t) + LHE (8)} C). 


pq +r. 


Dacă F(s) este improprie (m > n), partea polinomială nu mai este o simplă constantă 


(a cărei transformată Laplace inversă este proporțională cu 65), ci un polinom 


a cărui 


transformată Laplace inversă este o sumă de semnale singulare (impulsive). Această 


situație este în cadrată de Teoria Distribuțiilor și nu face obiectul acestui curs. 
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Exemplu. Să se calculeze funcția pondere a sistemului cu funcţia de transfer 


ERE e t3 


H -A 
(5) 2s2 + 6s +4 


Trebuie să determinăm £L7+(H(s)). Se observă că 


H(s) = 0.5+ ESE ee a 


Dao BIO a Gueta 


de unde 
h(t) = 0.55(t) + 0.5e-t —0.5e-%. 


Exercițiu. Demonstrați că L16(t))(s) = 1. 
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2. Stabilitate 


Ideea: un sistem este stabil dacă răspunde în mod adecvat unui stimul extern. 


STABILITATE: cerință fundamentală în proiectarea sistemelor de reglare automată 


(SRA). 


- SRA: ieșirea y “reproduce” referinţa r; forma de variaţie a lui y = forma de variaţie a 
lui r. 


- trebuie menținut controlul variaţiei mărimilor de ieșire; în plus, mărimile fizice care 
intervin nu pot fi oricât de mari, pentru că s-ar putea distruge instalaţia. 


- cerințele de proiectare impun realizarea unui regim staționar (permanent). Acest regim 
este realizat de către sistemele stabile. 
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Definiţia stabilităţii 
Reamintim că un semnal f : R — R este mărginit dacă există M > 0 astfel încât 


f(t)] <M, vteR. 


Un semnal măriginit are atât norma sup cât și norma œ finite. 


Definiţia 25. Un sistem y = T (u) este stabil în sens intrare mărginită- ieșire mărginită 
sau stabil BIBO, dacă pentru orice intrare u mărginită ieșirea rezultată, y=T(u), este 
de asemenea mărginită. 


Când este stabil un sistem de convoluţie ? 


Condiţii asupra funcţiei pondere și/sau asupra funcţiei de transfer. 


3Bounded Input/Bounded Output 
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Stabilitatea sistemelor de convoluţie 
Teorema 20. Fie X un sistem de convoluție y = h x u. Sunt adevărate următoarele 


afirmaţii: 


(i) Sistemul > este stabil dacă şi numai dacă funcţia pondere este absolut integrabilă 
pe R, h € LI(R), adică 


+00 
In. = RTA 


OO 


(îi) Dacă © este stabil atunci 


suply(t)| < |All suplu(?)|. 
teR teR 
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Stabilitate în sens nestrict 


Observaţia 27. Stabilitatea BIBO a sistemelor de convoluţie, caracterizată de condiţia 
h € LI(R), se mai numește și stabilitate externă în sens strict. 


Definiţia 28. Un sistem de convoluție y = h x u este stabil extern (sau BIBO) în sens 
nestrict dacă h(t) este mărginită. 


Dacă y = h x u este un sistem de convoluţie care admite funcţie de transfer rațională, se 
pune problema caracterizării stabilității (în sens strict sau nestrict) în termenii funcției 
de transfer a sistemului. 


Fie H(s) o raţională proprie, cu coeficienţi reali și ireductibilă, având polii p;, i = 1, k. 
In ce condiţii 


k 
h(t) = LHH (8) (t)= X m(t) ei + H(oo)ă(t) 
i=1 
este absolut integrabilă sau mărginită pe R ? 


Răspunsul: teorema următoare. 
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Stabilitatea sistemelor cu funcţii de transfer raţionale 


Teorema 29. Fie > un sistem având funcția de transfer H(s). Presupunem că H(s) 
este ireductibilă. Sunt adevărate următoarele afirmaţii: 


(i) Sistemul » este stabil BIBO în sens strict dacă şi numai dacă toţi polii funcţiei 
de transfer H(s) au partea reală strict negativă, adică 


Rep; < 0, Vi=1l,k. 


(ii) Sistemul X este stabil BIBO în sens nestrict dacă şi numai dacă toți polii 
funcției de transfer H(s) au partea reală negativă, adică 


Rep, < 0, T= l,k 


tar acei poli cu Rep; = 0 sunt poli simpli. 
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Exemple 


1. Fie H(s) = Polii sistemului sunt pı 2 = —2 +j“. Cum Rep <0, 
(s) PEEN] | | P1,2 > &J-3 P1,2 
rezultă că sistemul este stabil (în sens strict). 
| 25 — 1 no. . m 
2. Fie H(s) = Polii sistemului sunt pı 2 = +j și ps = —1. Deoarece 


(s2 +1)(s+1) 
Rep; < 0, și —j, j sunt poli simpli, rezultă că sistemul este stabil (în sens nestrict). 


Notă. In analiza și mai ales în sinteza sistemelor automate, suntem interesaţi de 
stabilitatea în sens strict. 


Terminologie: prin stabil înţelegem stabil în sens strict ! 


Cum testăm stabilitatea unui sistem oarecare ? In exemplele de mai sus am calculat 
ușor polii sistemului, iar în cazul general acest lucru nu mai este posibil. 


Condiţii asupra coeficienţilor. 
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Criteriul de stabilitate Hurwitz 


Problemă. Fiind dat un polinom oarecare p(s) cu coeficienți reali, în ce condiţii 
(exprimate cu ajutorul coeficienţilor) are polinomul toate rădăcinile în 


Co =iseC : Res <0}? 


Teorema 30. Fie polinomul p(s) = ans” +an—18”7! +. ..+a1s+a0, Cu an > 0. Atunci 
p(s) are toate rădăcinile în CT dacă şi numai dacă toţi minorii matricii Hurwitz (28) 
sunt strict pozitivi, adică 


ün- ün=3 Üp=E ars U 
| Ün ue a sxs Ü 
H”! > 0 vi=1: n, unde H := 0 Ün-i Qn-5 +i Epon (28) 


0 ün n=? 
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Observaţii 


Observaţia 31. 


1. Un polinom care are toate rădăcinile în C7 se mai numește și Hurwitzian (sau 
polinom Hurwitz). 


2. Ipoteza an > 0 nu restrânge în nici un fel generalitatea; dacă an < 0 se utilizează 
criteriul pentru —p(s) (care are aceleași rădăcini ca și p). 


3. Dacă există un coeficient nul a = 0 sau 2 coeficienți de semne contrare, aja; < O, 
atunci p(s) nu este Hurwitzian. 


4. Criteriul Hurwitz nu furnizează nici o informaţie cu privire la plasarea polilor - gradul 
relativ de stabilitate. 
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Exemplu: un polinom de gradul 3 


Fie p(s) = a3s% + a28? + as + ap, az > 0. Avem n = 3 și 


do QQ 0 
H := a3 di 0 
0 do QQ 


(1) HË! = az > 0. 
(2) HPI = ajaz — agao > 0. Cu (1) și (3) rezultă a, > 0. 


(3) HBI = ago(aia2 — agao) > 0. Cu (2) rezultă ao > 0. Aşadar a; > 0, i = 0,3 și 
a102 — a3a0 > 0 sunt condiții necesare și suficiente pentru ca p(s) să fie Hurwitzian. 
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3. Regim Permanent si Tranzitoriu 


3.1 Regimuri de Funcţionare Fundamentale. 


î.C PH 


Figura 21: Sistem de reglare automată 


Obiectivul reglării: procesul de a impune ca anumite variabile specificate (y) ale unui 
sistem să urmeze anumite evoluţii impuse (r), în prezenţa diferitelor perturbații (d). 
Altfel spus, ieșirea y trebuie să urmărească referința r. Trebuie “controlate” variațiile 
mărimii de ieșire. 


Forma de variaţie a mărimii de ieșire y permite aprecierea regimului de funcţionare al 
SRA, aceasta fiind direct influenţată de 

a) cauze externe (referință, perturbații) 

b) cauze interne (transformarea unor energii acumulate în sistem - ține de elementele 
constructive ale acestuia) 
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Regim forțat: forma de variație a lui y reproduce forma de variaţie a lui r (sau a lui d). 
Regimul forțat poate avea un caracter permanent sau staționar dacă mărimile exogene 
își păstrează forma de variaţie suficient de mult timp, caz în care el se va numi regim 
permanent sau staționar. 


Semnale cu caracter permanent sunt, de exemplu, semnalul de tip treaptă sau semnalul 
armonic. 


Atunci când o mărime exogenă (referință sau perturbaţie) își modifică forma de variație 
în timp (de exemplu, un semnal de tip treaptă este înlocuit cu un semnal armonic), 
atunci sistemul poate trece din regimul staționar de tip treaptă (1) într-un regim staționar 
de tip armonic (II). 


Regim tranzitoriu: regimul de funcționare al unui SRA aflat în trecerea de la un regim 
permanent la alt regim permanent. 


“Exemplu”:  TRANZITIA (de acum celebră) a societăţii românești de la un regim 
permanent (regimul Ceaușescu) la un alt regim permanent (?) 
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Sistemul nu poate urmări imediat schimbarea formei de variaţie în timp a mărimilor 
exogene: inerție datorată energiei acumulate în interiorul sistemului. 


Regimul liber: cauze interne (detalii complete - Teoria Sistemelor II). 


Pentru atingerea obiectivelor reglării trebuie realizat un regim permanent. 
Sistemele în care se poate realiza un astfel de regim sunt cele stabile. 


In această situație, componentele liberă respectiv tranzitorie (datorate unor cauze 
interne) se “sting” după suficient de mult timp. 


Instabilitatea face ca o parte din mărimile fizice care intervin în sistem să “o ia razna” 
sau să scape de sub control. 


O durată prea mare a regimului instabil poate conduce chiar la distrugerea 
instalaţiei / procesului. 


Morala: Stabilitatea - cerința fundamentală în proiectarea SRA. 
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Semnale cu caracter permanent 


Spunem că un semnal u(t) are caracter permanent dacă Jim ut 
— DO 


sau dacă lim u(t) nu există. 
t— o 
Semnale (standard) cu caracter permanent (sau persistente): 


t(n- 1) | 
ult) = 10), ult)= m I u(t) = Ae*. 


Observaţia 32. Polii transformatelor Laplace ale unor astfel de semnale sunt în 
(se C : Res > 0}. 


De exemplu, u(t) = t1 (t) este un semnal persistent (lim u(t) = +oo), cu L{u(t) ys) = 


1 
-z pe când v(t) = e”! 1(t) satisface lim v(t) = 0, cu L{v(t)} (s) = , Și nu are 
s => oo 


caracter permanent (se “stinge” după suficient de mult timp). 


s+ 1 
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Răspuns permanent și răspuns tranzitoriu 


Fie sistemul STABIL 


având polii p; cu Rep; <0,2=1l:k. 


Considerăm semnalul exogen (intrare, referință, perturbație) e(t) a cărui transformată 


rels) 
Pe(s) 


Hels) pej, verifică Repe; > 0, J=l:n. 


Laplace este He(s) = „ Presupunem că e(t) are caracter permanent, deci polii lui 


Rezultă 
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Cum p și pe nu au rădăcini comune, Y (s) se descompune în fracţii simple astfel“: 


Ya = SO A 


Deducem că E 
y(t) = 5 A, ePit + ` Aej ePeit (29) 
i=1 j=1 
yt(t) Yp(t) 


Componenta y+ a răspunsului se numește răspuns (regim) tranzitoriu iar componenta yp 
a răspunsului se numește răspuns (regim) permanent. 


Se observă imediat că Jim y+(t) = 0 și că yp(t) reproduce forma de variație a semnalului 
—> OO 


exogen e(t), având un caracter permanent (explicaţi de ce). 


4Pt, a fixa ideile, am presupus că atât H(s) cât și He(s) au numai poli simpli. Așa cum se știe, dacă apar multiplicităţi, 
în locul coeficienţilor A;, Aej apar în expresia (29) polinoamele r;(t), rej(t). 
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3.2 Răspunsul Sistemelor la Semnale de Intrare Standard 


1. Intrare de tip treaptă: u(t) = 1(t); He(s)=-—. 


Rezultă n 
t) = A; ePit + H0) 1(6). 
u(t) = Š 4i eP + HO) Lt) (30) 
= Yp(t) 
yt(t) 


Răspunsul la un semnal de intrare de tip treaptă se mai numește și răspuns indicial. 


Observăm că (/(s) stabilă și strict proprie) 


(00) lim sY (s) = H(0), y(0) "lim sY (s) = H(œ0) = 0. 


S —> OO 
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pnl 


2. Intrare de tip polinomial: u(t) = EN 1(t); Hels) = —. 


unde 


Rezultă 


s E (0) Ea, 
=l 
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Ss — jw 
k A, Ao 
Y(s) = H(s i Ao = sY (s) s=; = H (j C 
(6) = H6) = tag 5Y Hu) e 
Rezultă r 
[= A, ePit + H(j Jor, 2 
J= SAP + HOw) e (32) 
Ti Yp(t) 


yt(t) 


Așadar, răspunsul permanent sau staționar al unui sistem stabil la o intrare de tip 
armonic este un semnal armonic de aceeași frecvență cu semnalul de intrare, având însă 
amplitudinea și faza modificate. 


Reamintim că H(ju) = |H (jw)| ei sl") este răspunsul în frecvență al sistemului. 


Mai multe detalii: răspunsul sistemelor elementare la intrări standard. 
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3.3 Performanţe de Regim Permanent și Tranzitoriu 
Revenim la schema standard de reglare din figura 21. 


1. Cerinţe relative la regimul permanent sau staționar. 
e Regim permanent la semnal de tip treaptă unitară. 


IDEAL: eroare staționară ZERO, (oo) = 0. In acest caz, y(00) = r(oo) = 1. 
In practică, se va impune o eroare staționară foarte mică. In acest caz, y(0o0) ~% r(oo) = 1. 


e Regim permanent la semnal de tip rampă. 
Eroare staționară e(00) = r(o0) — y(o0o) suficient de mică. 
e Regim permanent la semnal de tip armonic. 


Amplificarea /atenuarea |H (jw)| și defazajul arg|H(jw)| nu trebuie să depășească anu- 
mite valori impuse, în intervale de frecvență specificate |u_ , w4]. 


Performanţele de regim staționar sunt specificate pt. cele 3 tipuri (uzuale) de semnale 
standard. Performanţele de regim tranzitoriu vor fi specificate doar pt. răspunsul la 
intrare de tip treaptă unitară. 
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2. Cerinţe relative la regimul tranzitoriu. In figura de mai jos este reprezentat grafic 
răspunsul la treaptă ya(t) al unui sistem stabil H(s). Reamintim că ya = yt + Yp, unde 


Yp(t) = H(0) 146). 


Va 
“max þ----------------,5 


1.02 H(0) 


HOLA] i 
0.98 HO- Sa > 
0.9 HO) 


Figura 22: Regim tranzitoriu la intrare treaptă unitară 
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e timp de răspuns suficient de scurt 


- timp de creștere te: timpul în care ya(t) crește de la 0.1 H(0) la 0.9 H(0) (timpul 
necesar ajungerii în vecinătatea unui nou regim staționar). 
- timp de vârf t,: timpul necesar pt. ca yı(t) să atingă valoarea maximă. 


e calitatea urmăririi treptei 


- timp tranzitoriu tę: timpul necesar încadrării răspunsului yı(t) în intervalul 
[0.98 H(0) , 1.02 H(0)| (timpul necesar componentei tranzitorii să devină aproape 
zero) - marchează terminarea regimului tranzitoriu. 


Așadar |y:(t)| < 0.02 H(0), pentru t > te. 


U1,max a Yalt) — Y1,max as H (0) y 


- suprareglajul o: o = 
Yp(t) H (0) 


- valoarea de vârf: Y1,max- 


Studiu și calcul detaliat: sisteme de ordinul II. 
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3.4 Răspunsul Sistemelor Elementare la Intrări Standard. 


Considerăm sisteme de ordinul |, respectiv de ordinul II. 


a) Răspunsurile sistemului de ordinul | la intrare de tip treaptă, rampă și de tip armonic. 
Fie 


K 
H(s) = ——, 
(5) T's+1 
unde K,T > Q sunt factorul de amplificare respectiv constanta de timp a sistemului. 
. Z K K K J 
- Răspuns la intrare treaptă: Y (s) = ———— = — — r- Rezultă 


s([3+1) S Saa 


-1 
v(t) = (K — Ket) 10). 
Regim staționar: identic cu intrarea (pt. K = 1). 


Suprareglaj: o = 0. 
Timp tranzitoriu: t; ~ 4T (|=Ke-7| < 0.02K, Vt > 4T): 
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Figura 23: Răspunsul la treaptă al sistemului de ordin | 
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l K KT H(0)  H'(0 
- Răspuns la intrare rampă: Y(s) = 27531) DTF L ) T A 
s? (Ts s+ S S 


Yramp(t) = KT e77 1(t) + K(t — T) 1(¢); 


Regim staționar: NU este identic cu intrarea (pt. K = 1); are doar aceeași formă de 
variație. 
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Figura 24: Răspunsul la rampă al sistemului de ordin | 
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| K 
- Răspuns la intrare de tip armonic u(t) = et: Y(s) = = 
P Ni g (5) (Ts + 1)(s — jw) 
H(j — H(} i 
CT (jw) I 
K | 
ya = T e™ t>o, 


T k P AR e T + 57 a E 

Lpa lLpjwr 
Regim staționar: Semnal de aceeași frecvență, cu amplitudine și fază modificate. 
Rey~(t) - răspunsul la intrare coswt1 (t) 


Imy~n(t) - răspunsul la intrare sin wt1 (t). 
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b) Răspunsul sistemului de ordinul II la intrare de tip treaptă. Fie 


82 + 2Cums +2! 


H (5) 


unde wn > 0 și 0 < Ç < 1 sunt pulsația naturală, respectiv factorul de amortizare. 


Polii sistemului H(s) sunt complex conjugați: 
t | 
pp => id E] Way p= (7 = Oq E ] ua. 


- Răspuns la intrare treaptă u(t) = 1 (t) : 


Wn H (0) S— Oq Wd a 


Y = — v — 
(5) s(32 + 2Cups + w2) S (s= 04)? twa (5-04)? + wi Vi- e 


Deducem că 


yı(t)=1 Aa sin(wn y 1 —(2t+gp), t>0; cosy =¢ 
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Figura 25: Răspunsul sistemului de ordinul Il la intrare treaptă unitară 
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G = 0: y(t) = 1 — cos wnt - regim oscilant. 
(= 1: y(t) = 1 — wpte” ont — e-ent - regim critic. 
C > 1: regim supracritic - sistemul are 2 poli reali simpli. 


C < 1: regim subcritic - cazul care ne interesează; apar oscilațiile amortizate. 
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Performanţe de regim staționar 


Determinăm punctele critice ale funcţiei y1(t): 


d t n . ae 
y(t) ae sin(wn y 1 — °t) e Swnt — Q, 


dt vi- 2 


Rezultă ty = ——=——, de unde 


| —e-Sntk k= par minime locale 
_ 1 _ (_1ko—Cuntk — p 
Wid = =e | l+e-ntk k= impar maxime locale 


_ ST 
Se observă că ymax = y(ti) =1l+e V1, deci suprareglajul este dat de 


_ 6T 


o=e VË, (33) 
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Notă: Suprareglajul depinde doar de factorul de amortizare”. 


o% 


100% 


1 


Figura 26: Dependența suprareglajului de factorul de amortizare 


SErată pe grafic: pt. C = 0.4 suprareglajul este cca. 25%-26% 
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îi ui 0.2 04 05 0.6 0.7 
Tabel cu valori estimative o 50% 25% 16% 10% 5% 


Estimarea timpului tranzitoriu t+. 
Dacă pt. un extrem local |y+(tk)| < 0.02, atunci |y,(t)] < 0.02 pentru orice t > tz. 
Se verifică (la fel ca la estimarea t; la sistemele de ordin I) inegalitatea |y:(t)| < 0.02 


. 4 < 
pentru orice tg > —. Rezultă 
n 


4 
lyelt)| < 0.02, vt > ta. Rd ty: 


n 


Pentru 0.3 < C < 0.8, timpul de creștere poate fi apreciat cu ajutorul formulei 


om A 3 
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Timpul de vârf, tẹ, este chiar ti = 


Plasarea polilor în raport cu cerinţele de performanţă ale regimului 
tranzitoriu 


Figura 27: Plasarea polilor sistemului de ordin Il 
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Problemă: Se cere determinarea zonei din planul complex în care trebuie plasați polii 
unui sistem de ordin Il, pı 2 = —Çwn £ j wny/ 1 — G?, astfel încât regimul tranzitoriu să 
în deplinescă anumite cerinţe de performanţă impuse: o < oo, tt < tto, te < teo. 


e Suprareglaj. o < oo implică C > (o (vezi figura 26); cum ( = cosp (0 < p < 1/2), 
rezultă p < po, po = arccosco. 


Po 
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e [imp tranzitoriu. Avem t; = a < to, de unde Cwn > =- 
n tO 
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| ` 1.8 
e Timp de creștere. La fel ca mai sus, se deduce că wp > —. 


cO 
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Toate constrângerile 
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Reglarea poziţiei unghiulare a unei antene de satelit 
Exemplul 33. Ecuația mișcării în plan vertical a unei antene de satelit este dată de 


JÖ +b = M 


unde 0 este unghiul de înălțare a antenei, J = 6. 10” kg. m? momentul de inerție al 
acesteia, b = 2 - 104 N. m - sec coeficientul de frecare iar M cuplul motor. 

Se dorește aducerea antenei intr-o poziție (unghiulară) specificată, 0,, utilizând o schemă 
de reglare în reacție inversă cu un compensator de tip proporţional (K), astfel încât 
suprareglajul SRA să nu depășească 10% iar timpul de creștere să nu fie mai mare de 
&Usec. 


Figura 28: Sistem de reglare automată a poziției unghiulare a unei antene 
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Notaţie: u = M (comanda), y = 0 (ieșirea), r = 0, (referința). 
Verificaţi că funcţiile de transfer ale sistemului, respectiv ale sistemului în buclă închisă, 
sunt date de 


1 K KJJ 
(5) Js? + bs’ (5) Js +bs+K s2+b/Js+K/J 


|K b 
Sistemul în buclă închisă = sistem de ordinul Il, cu wn = 4/ — și C = ——. 
J IYIK 


Din graficul 26 și tabelul aferent se deduce că o < 10% implică ¢ > 0.6, deci 
VK < b/(1.2VJ). Rezultă K < 700. 


Deoarece te < 80sec, se obține imediat inegalitatea w, > 1.8/80 și rezultă 
VK > VJ1.8/80, de unde K > 300. 


Concluzie: Pentru 300 < K < 700, sistemul de reglare este stabil, are eroare staționară 
zero, suprareglajul mai mic de 10% și un timp de creștere care nu depășește 80sec. 
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4. Reprezentarea in frecventa a sistemelor 


Electronică: curba de răspuns în frecvență a unui amplificator arată modul cum se 
schimbă câștigul unui amplificator atunci când frecvența semnalului de intrare variază. 


Banda amplificatorului: plaja de frecvenţe unde amplificarea este practic constantă, la 
un nivel apropiat de valoarea sa maximă. 


Răspunsul în frecvență al unui sistem care admite funcţie de transfer H(s) este 


H (jw) = H(5)ls=io 


Pe de altă parte, H(jw)e'“! este răspunsul permanent al unui sistem stabil la o intrare de 
tip armonic, u(t) = e“. Mai precis, dacă la intrarea sistemului avem u(t) = Asin wt, 
răspunsul permanent (sau staționar) al acestuia va fi 


Yp(t) = AJH (jw)|sin(wt + ¢w)), 
unde ġ(w) = arg| H (jw). 
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Reprezentări în frecvenţă 


Se numește reprezentare în frecvență a funcţiei de transfer H(s), o reprezentare grafică 
a dependenței numărului complex H(juw) de w, pentru w E€ R. 


Numărul complex H(jw) poate fi reprezentat fie în forma 
i à : š not : 
H (jw) = Re(H (jw)) + jlm(H(ju)) = U(w) + JV (w) (34) 
fie în formă exponențială (sau polară) 


H(jw) = |H (jw)| ial E Hu) ei (35) 


Vom introduce 2 tipuri de reprezentări: 
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L Caracteristicile amplitudine-pulsație și fazā-pulsație sau caracteristicile 
(semi)logaritmice sunt graficele amplitudinii H (w) și respectiv ale fazei (w). Aceste 
caracteristici se mai numesc și diagrame (de tip) Bode. 


2. Locul de transfer sau locul Nyquist este hodograful lui r(w) = (U (w), V (w)) 
pentru w € R, adică submulțimea lui R? = C definită de {(U (w), V (w)) |w € R}. 


Există și alte reprezentări, cum ar fi caracteristica amplitudine-fază sau locul de transfer 
invers. Aceste caracteristici se definesc în mod similar și pentru sistemele discrete (detalii 
în semestrul Il...) 


Observaţie. Aceste reprezentări pot fi obținute pe cale experimentală, fără ca funcţia 
de transfer a sistemului să fie cunoscută! 

Există instrumente care măsoară variațiile de amplitudine și fază în în raport cu 
modificările frecvenței semnalului sinusoidal de intrare. 


Se poate deduce chiar modelul unui sistem utilizând astfel de dispozitive. 
Hewlett-Packard 3562A Dynamic Signal Analyzer. 
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Deoarece H(s) este o rațională cu coeficienţi reali, rezultă că 


H(ju) = Ul) — jV (w) = Hw) = H(-ju) = U(-u) + jV(-u) 


de unde U (w) = U(—uw), -V (w) = V(—uw) pentru orice w € R. Cu alte cuvinte, U este 
o funcţie pară iar V o funcţie impară de w. Cum 


H(w) = VU? (w) +V? (w) și Șlu) = arcte 


rezultă că H este o funcție pară iar g o funcție impară de w. 


Așadar este suficientă trasarea caracteristicilor pentru valori pozitive ale pulsației w. 
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Diagrame de tip Bode: H(u), (u) 


Aceste caracteristici se mai numesc și (semi)logaritmice deoarece pe abscisă se alege 
o scară logaritmică pentru w, x = loguw, iar pe ordonată y = H(u) și y = p(w) se 
măsoară în decibeli 


[Ha := 20log H, H >Q, 
respectiv în radiani (sau grade). 


A IYlag 


i decadă i 
> 


L L > 
0 101 102 logx 


Figura 29: Abscisă logaritmică: o decadă 


log 2 ~ 0.3, log5x:0.7, log 3 2: 0.48 (se poate considera chiar 0.5) 
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Reprezentarea sistemelor elementare 


K 
1. Elementul (multi)integrator: H (s) = n K >0,gqeZ. 
S 


Aşadar [H(w)ļaB = LA Jag ace 20q log. w. 
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Bode Diagram 


Magnitude (dB) 
DD 
O 


90 


Phase (deg) 
O 
T 
| 


107 10! 10° 10" 10° 


Frequency (rad/sec) 


Figura 30: Multiintegrator: q = —2, —1, 1,2 
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1 
2. Element de ordinul 1: H (s) = TIT T >O: 
S 


1 1 


Hiu =. H = 
Go) = pp H4w) TEN 


olw) = arctg(—wT). 


Rezultă că 


[H (w)l]as = 20log (1 + b2T2)-î = —10 log (1 + w2T2). 


1 22 
a) joasă frecvență: w << T' 1 +w Ti sl. 
In acest caz, avem |H (w)lag ~ 0, e(0) = arctg(0_) = 0. 
l 22 272 
b) înaltă frecvență: w >> T’ 1 +w Ti x w TA. 
In acest caz, avem [H (w)]ag ~ 20log 1/7 — 20 logw, (+œ) = arctg(—œ) = —3. 
1 


c) medie frecvență: w = r 
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In acest caz, avem [H (w)]ag = —10log2 = —3dB, (1/1) = arctg(—1) = —3. 


1 l r 
wp = T se numește pulsația de frângere. 


Bode Diagram 


Magnitude (dB) 
o 


Phase (deg) 
A 
a 


_90L PI oi biarin i bi == 
10 107 10° 10 10 
Frequency (rad/sec) 


Figura 31: Element de ordinul 1 
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w2 


i | o E 1 
2. Element de ordinul 2: H (s) = CEDETTE = CTR FL 0 si 4 < iA 
w? 1 
A E = NN e ee NI 
(d) —w? + 2(unu +w}?  1—x2+j 2 Uz) 
unde z = 5- Ro 0) pulsație normată. Avem 
232 2.213 —2çx 
(0) = HGE) = [PP +4], go) = arcte s 


Rezultă că 


LH (wlan = —10 log (1 — 22)? + 4622 


a) joasă frecvență: w << wn, x << 1. (1-2): + 40272 x 1. 
In acest caz, avem |H(w)lap x —10log 1 = 0, (0) = arctg(0_) = 0. 
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b) înaltă frecvenţă: w >> op, 2>> 1. 1— g?) +4 r? x rt. 


In acest caz, avem |H(w)ļlap œ~ —40logx = 40logwn — 40loguw,  g(+oo) = 
arctg(0_) = —r. 


c) medie frecvență: w =uw,, z=1. 
In acest caz, avem |H(w)lag = —10log4€? = 20 log ze =: [ze las (un) = 


arctg(—00) = —3. 


Puncte de extrem. Determinăm punctele critice ale funcţiei |H(w)laB(z): 


ALA (w)]an(z) 


=0 & za =0saul—2€2—a2=0. 
dr 


Rezultă z — 0 și, dacă € < 1/y2 — 0.707, x — 1 — 262. Se poate arăta că acestea 
sunt puncte de maxim, valoarea amplitudinii (exprimată în dB) fiind O în primul caz și 
respectiv 

[Hasau lag = —10 log 4¢°(1 E a 


în cel de-al doilea. 
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Bode Diagram 


20 |- 


Magnitude (dB) 


-90 F — 


Phase (deg) 


-135 | =] 


-1806n s us al , ni E 


Frequency (rad/sec) 


Figura 32: Elemente de ordinul 2 cu w, = 3; parametrizare după Ç. 
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Concluzii procedurale 


Se trasează mai întâi caracteristicile asimptotice: 


0 0O<u<z 


Sistem de ordinul 1: LH A(w)las = l —90 anu? i sy 
T 


Sistem de ordinul 2: [H4 (w)ļaB = | DN 


Se deduc valorile caracteristicii exacte la pulsațiile de frângere 1/T, respectiv wn, și se 
trasează caracteristica exactă (ținând cont de valorile lui Ç pentru sistemele de ordin 2). 


De exemplu, dacă ( = 0.5 < 0.707, caracteristica amplitudine-pulsație admite un maxim, 
iar [H (wn)]ag = —10log 4¢? = 20 log ze = [zzozlaB = 0. 


a : l T l 
Variația totală a fazei este de E pentru elementul de ordinul 1, respectiv de —r pentru 
cel de ordinul 2. 
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Trasarea calitativă a caracteristicilor de tip Bode 


Fie funcţia de transfer (strict) proprie 


1 geZ, KER, 7(0)=35(0)=1. (36) 


r(s) _ K 7(s) 
s1 p(s) 


q se numește tipul funcției de transfer iar K = H (0) este factorul de amplificare/atenuare 
sau câștigul (DC gain). 


Punem în evidență factorii specifici care apar în expresiile lui 7(s) și p(s): 


Presupunem (deocamdată) că toate constantele de timp sunt pozitive (Tę > 0), că 
0O<(<lș A >0. 
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In expresia lui H(s) apar 4 tipuri de factori: 


Factorul constant K. Caracteristica amplitudine-pulsaţie este o constantă 
iar faza este egală cu 0. 


Poli/zerouri în origine (jw). 
Poli/zerouri reale (jwT + 1). 


Poli/zerouri complex conjugate (1 — w° T? + j2(Tuw). 
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Notă: Dacă H(s) = H,(s)Ho(s), atunci 
|H (w)]as = [H1 (o)lag + |Ha(o)lag 
arg|H (jw)] = arg Hi (jw)] + are] Hu). 


In consecință, 
a n n! 
(H(o)la = = + UL + juTrilag +X C1- 0273, + j2GiTrivla + 
dB i=1 i=1 


n, n”! 
-X _ |1 + jwfpjlas — X1 — 272; + 326; TpiwlaB 
jJ=1 j=1 


și caracteristica amplitudine-pulsaţie se obține adunând caracteristicile fiecărui factor în 
parte. 
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Similar, deoarece 


/ AA 


Ny Ny 
T a | 
olw) = E + 5 arg(1 + jwTri) + 5 arg(1 — ww2T2 + j2C; Tri) + 
i=1 i=1 
ni, T 
-X arg(1 + jwTp;) — > arg(1 — w Th + j2; pi), 
jJ=1 jel 


obținem caracteristica fază-pulsație însumând caracteristicile fiecărui factor elementar în 
parte. 
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Procedura de trasare 


1. Se pune sistemul în forma (36). 


2. Se determină valoarea lui K in dB, [K]ap = 20log K și se marchează punctul 
(1, KlaB). 


îi Í 
3. Se așează în ordine crescătoare, pe axa logw, toate pulsațiile de frângere węg = T 
s YA s i A k 

Se trasează câte o verticală în fiecare pulsație de frângere. 


4. Se trasează asimptota de joasă frecvență trecând prin punctul (1, |K]ag) și având 
panta de —20q dB/decadă. 


5. La intersecția caracteristicii asimptotice cu fiecare verticală, se modifică panta acesteia 
cu +20 sau +40 dB/decadă, după cum respectiva pulsaţie de frângere corespunde unui 
factor de ordinul 1 sau de ordinul 2, plasat la numitor sau numărător. 


6. Caracteristica fază-pulsaţie se trasează prin însumarea caracteristicilor fiecărui factor 
elementar. 
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Exemple 


5 1+0.1ls 


H(s) = 3 F058) (14 0.65/50 + (3/50) 


K=5-q= l; [las = 20log 5  14dB. 


Ti = 0.5, W1 = 2: To = 0.1, W2 = 10; 13 = 0.02, W3 = 50. 


Se trasează caracteristica asimptotică în conformitate cu regulile expuse la 4 și 5. 
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40r 


wi w2 w3 


-100 1 o i | a Ea i a 
10 10 10 10 10 


Figura 33: Caracteristica asimptotică 
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Caracteristicile “reale” 


Bode Diagram 
50 m = 


Magnitude (dB) 
| 
(S 
oO 


-100 


-270 E , r a i i p a i giii i i iii 7 
107 10 10' 10 10 


Frequency (rad/sec) 


Figura 34: Caracteristicile amplitudine-pulsaţie și fază-pulsaţie 
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Alte situaţii 


1. K <0. 


- H(w) nu se modifică 
- la expresia lui (w) se adună arg| K| = —r (vom conveni că arg|—1] = —7). 


2. T < 0 (element de ordinul |). 


j E 1 LoT, 
1+juT  1—jo|T| 1+w?T?’ 


- H(w) nu se modifică (explicaţi de ce) 
- Pl) = arctg(uo|T]) = —arctg(—w|T]) 


Concluzie: 
arg[(1 + juT)"1] = sgn T arg[(1 + ju|T])) 
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3. 1 <0,0< Å< 1 (element de ordinul II). 


| 1 tea 070 e 
PT <O, HG2) = a Pofta 0 a + ară 


- H(w) nu se modifică (explicaţi de ce) 


- ġlw) = arctg 2/9 = —arctg í ele) 


] — z2 | 2 


Concluzie: 


arg](1 — 272 + j26wT)1] = sgn T arg|(1 — 272 + 2cw|7])] 


Notă: situația T > 0, —1l < Ç < 0 se tratează în aceeași manieră. 
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Sisteme de fază minimă 


Fie 
o 10 1 — Tis 


10 1+ Tis 
i i e IO EI a a nn l 
(s) TETT 


a ale 
9) = ea e, 
unde 1, > 15 >Q. 


Ne interesează caracteristica fază-pulsaţie a celor 2 sisteme (vezi figura 35). 


Defazaj total H(s): —3r/2. 
Defazaj total G(s): —r/2. 


Concluzie: Sistemele cu zerouri nestabile (care se află în semiplanul drept, Re z > 0) au 
un defazaj total mai mare decât sistemele cu zerouri stabile (Rez < 0). 


Sistemele (stabile) cu zerouri stabile = sisteme de fază minimă. 


Sistemele de fază neminimă sunt greu de “controlat” (comandat). 
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Bode Diagram 
100 T Tr] T FT EE E E T TIT TITI T 


50 


Magnitude (dB) 


—100 


-158 


Faza lui G(s) 


] 
[te] 
O 


Faza lui H(s) 


Phase (deg) 


-180 - 


970 ha atei pje fe jaf iu, i 
10° 10° 10! 10% 10 10 10% 


Frequency (rad/sec) 


Figura 35: Sistem de fază minimă 
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Locul de transfer (locul Nyquist) sau hodograful 


Reamintim că locul de transfer sau locul Nyquist al sistemului H(s) este hodograful 
lui T() = (U(w), V (w)) pentru w € R, adică submulțimea lui R? = C definită de 
((U(w), V(u)) |w e RI, unde 


U(w) := Re(H(jw)); V(u) = Im(H(jw)). 


De asemenea, am arătat că U (w) este o funcție pară de w iar V (w) este o funcție impară 
de w. Este deci suficient să trasăm locul de transfer pt. w € [0, c0). 


Observația 34. (importantă !) 

Hodograful este reprezentarea conformă a axei imaginare s = jw din planul complex 
al variabilei s, definită de funcția analitică H (s). Axa imaginară se reprezintă ca un 
contur închis (conturul Nyquist), capetele de la infinit fiind legate printr-un semicerc de 
rază R — 00. Așadar, hodograful va fi la rându-i un contur închis în planul complex 
H (s). 
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Locul de transfer al sistemelor elementare 


K 
Integratorul : H(s) = —, K > 0. 
S 


Tabelul de variație: 
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Nyquist Diagram 
10 T T T T 


Imaginary Axis 
e 
| 
| 


-1 -0.8 —0.6 -0.4 -0.2 0 0.2 0.4 
Real Axis 


Figura 36: Integrator - loc Nyquist 
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1 
Elementul de ordin |: H (s) = TLT LN 
S 


1 l= 000 1 wT 
H siu H k = a E 3, V = — 


w| o 1/T 
Tabelul de variație: [U| 1 12 Q| 


vo = o 


o. ae ză 2 _ 
Se poate elimina w din expresiile lui U și V: w272 = Er. Rezultă 


(EE E). 


care este ecuația unui cerc de centru (1/2,0) și rază 1/2. Se observă că hodograful 
sistemului de ordin 1 nu depinde de T 
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2 


“n Wn > 0,G € [0,1). 


Elementul de ordin Il: H(s) = = ————— 


2 
Wn 


R E E a FER 


_ z=w/wn 1 


aena H(jæ) = U(z) + jV(z); 


unde 


| — 2 24) 
(1 — z?) + 4622) E 


(a — g2)? + 42x2 


Atunci când Ç creşte (către 1), intersecția cu axa imaginară se apropie de valoarea —1/2. 
Când Ç scade (către 0), intersecția cu axa imaginară se apropie de —oo. 
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Reamintim faptul că pentru 0 < Ç < a modulul lui H(jw) prezintă un maxim la 
— v 1 . : 
frecvența w = wny/ 1 — 262. Dacă 5 < Ç< 1, atunci |H(jw)| scade monoton cu w. 


Imaginary Axis 
O 
Imaginary Axis 


L L L L ji L | L 
-1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 1.5 2 
Real Axis Real Axis 


(a) (b) 


Figura 37: (a) Elementul de ordin |; (b) elementul de ordin 2 
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Trasarea calitativă 
Comportarea hodografului la joasă frecvență: w — 0. 
Scriem sistemul în forma standard 


Ma Ri) unde 7(0) = p(0) = 
ii Aare sa pls) de F(0) = p(0)= 1, q EZ. 


H (0) = 0 dacă q < 1, H(0) = K dacă q = 0. 


Fie q > 1. Dezvoltăm în serie în jurul lui O: 


0 “i . — 
H(s) ~ K |S +a]; ao = 1. (37) 
Pt. q = 1 egalitatea (37) se scrie (s = jw) 
K K K 
H(jw) x a = Kay—=j 5, Uw = Ka, =const., V(u) = Îi: API ia 
ju w Ww w 
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Pt. q = 2 egalitatea (37) se scrie (s = jw) 


a a Ka , Ka 
H (jw) ~ Kat gu + aa — = + Ka+j (-==) = 
U(w) va 


Se poate deduce că 


care reprezintă ecuația unei parabole în planul (U, V). 


Acest tip de analiză se poate efectua și pentru q = 3. 
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v 
) 
IV 
A 
EE 
Y 
| 
Figura 38: Comportarea locului de transfer la joasă frecvență 
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Comportarea hodografului la înaltă frecvență: w — +oo. 
Fie 
(5) r(s)  bmsm+...+bis+ bo 
5) E = 

p(s) ans” +... +as+ ao 
și e := Olp(s)] — O|r(s)] excesul poli-zerouri al lui H (s). 


Pentru s — œ avem H(s) x Pm Ł, de unde 
w 1l 
Hua =- = — 0. 
Ga ju 


Rezultă i 
_ ze, 0 
p(00) o | T bm < 0 
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Tabelul de variaţie. 


Pentru o trasare cât mai precisă a locului Nyquist se realizează un tabel de variaţie a 
funcțiilor U (w) și V(w), evidenţiindu-se 


- intersecțiile cu axele (V (w) = 0, U (w) = 0); 
- posibile discontinuități, eventual asimptote oblice; 


- domeniile de valori ale funcţiilor U(w) = 0, V(w) = 0 precum și proprietăți de 
monotonie ale acestora. 
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Criteriul Nyquist 
r(5) 


Se consideră un sistem H(s) = —— conectat în reacție inversă: 


p(s) 


Figura 39: Sistem în buclă închisă 
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H 
Funcţia de transfer în buclă închisă: T'(s) = a 
S 


Problemă: Se cer condiții (necesare și suficiente) astfel încât sistemul în buclă închisă 
să fie stabil. 


Soluţie: Locul de transfer al lui H(s). 
Ideile: 
- polii sistemului în buclă închisă = zerourile lui 1 + H(s). 
- buclă închisă stabilă <> 1+ H(s) poate avea doar poli în Res > 0. 
- teorema variației argumentului (contur Nyquist & loc Nyquist). 
- polii lui 1+ H(s) = polii lui H(s). 


Rezultatul: 
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Teorema 35 (Criteriul Nyquist). 


In ipoteza că H(s) este ireductibilă (nu au loc simplificări poli-zerouri), sistemul în 
buclă închisă. este stabil dacă și numai dacă hodograful lui H(s) nu trece prin punctul 
critic —1 + j 0 sar variaţia argumentului vectorului cu originea în punctul critic şi 
vârful pe ramurile continue ale locului de transfer (acest vector este chiar 1+ H(s)!), 
atunci când w variază de la 0 la +oo (exceptând punctele de discontinuitate), este 
egală cu (2N, + 10) 


Aarg|l + H(s)lueto,o tu; | p(ju)=0) = (Np + qo)3. 


Aici N, este numărul polilor lui H(s) din semiplanul drept (Res > 0) iar qo este 
numărul polilor lui H(s) de pe aza imaginară (Res = 0), incluzând şi ordinele de 
multiplicatate. 
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Exemplu 


l 


Fie sistemul P(s) = 10633 


Să se analizeze stabilitatea sistemului de reglare de mai jos în funcție de parametrul K. 
r — e u 
[N 
Cp 


Figura 40: Sistem în reacție inversă 


y 


| K 
In figura 41, C(s) = —. 
S 
Funcția de transfer în buclă deschisă (sau pe calea directă) este ireductibilă, 


K 


H(s) = P(s)C (s) = GIGI 
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Trasăm hodograful lui H(s). 


Tabel de variaţie. Intersecţii le cu axele. 


Distingem 2 situaţii: când intersecţia cu axa U este 


P MITE: K 
- la stânga punctului critic: E: <=] s ULK 


- la dreapta punctului critic: = > —l 4 0L K. 
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K<6: sistem stabil K>6:sistem instabil 


Nyquist Diagram Nyquist Diagram 
1 T T T T 3 F T T T T T = 


0.8 | J 


E K2] 
x< x< 
< < 
2 > 
Ss Ss 
= £ 
D D 
S S 
£ £ 

-24 A 

-0.8 F — 
=] l j! i | l | l | Í 
-1 -0.8 —0.6 —0.4 —0.2 9) -3 -2.5 -2 -1.5 -1 -0.5 9) 
Real Axis Real Axis 


Figura 41: Sistem în reacție inversă 
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In figura 41 (a) variația argumentului vectorului 1 + H(juw), când w parcurge intervalul 
(0, 00), este 


m T 
A = arg|l + H(jw)lu=+o — arg[1 + H (jw)ju=0, = 0 — =. TJ 
In cea de-a doua situație (vezi hodograful (b) din figura 41) avem 
l . T ƏT 
A =arg|l + H(jw)lu=+%0 — argll + H(jw)lu=0 = (27) - (5) == 


Sistemul are N, = 0 (nici un pol în Res > 0) și qo = 1 (un pol pe axa imaginară, chiar 


în 0). Rezultă 
TOT 


2N 
( p + do) 2’ 


deci sistemul este stabil în buclă închisă pentru 


0<K<6. 
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Capitolul 4: CONCEPTE FUNDAMENTALE ALE BUCLEI DE 
REACTIE (FEEDBACK) 


Ce Dorim ? 

Bucla de Reactie 

Stabilitatea Interna a Buclei de Reactie 
Stabilizare 

Performante Asimptotice 


Solutia Problemei Reglarii 
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1. Ce Dorim ? 


Acest capitol este cel mai important pentru intelegerea scopului si mijloacelor auto- 
maticii ! (pica sigur la examen !!!) 


Introducem notiuni de baza cum ar fi conexiunea in bucla de reactie (feedback) si 
caracterizam principalele ei proprietati ce o transforma in instrumentul fundamental 
al automaticii (asigura stabilizarea + reglarea)! 


Punctul de plecare: Avem un sistem care nu face ceea ce dorim ! Cum procedam? 


Doua variante posibile: 

e || schimbam cu altul care face ceea ce dorim (dar oricum trebuie sa stim CUM sa 
cerem / construim sistemul nou...) 

e || modificam (CUM ?) pe cel pe care il avem deja. 


Analiza sistemelor din primele 10 Cursuri a incercat sa va familiarizeze cu diverse clase 
de modele sistemice si cerinte naturale asupra lor. 
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Clasa sistemelor considerate: 

e LINIARE 

e INVARIANTE IN TIMP 

e FINIT DIMENSIONALE 

e CAUZALE 

e O INTRARE — O IESIRE (SISO) 

Punctul de vedere: 

Intrare-lesire (1/0) (modelele sunt fctii rationale proprii) — 
Teoria matematica de baza: 

Functii de o variabila complexa 


Ce Dorim de la un astfel de Sistem ? 


Multe ... 

De exemplu: Sa fie stabil (semnalele externe de energie finita sa rezulte in semnale 
interne tot de energie finita), sa urmareasca clase de semnale prescrise in conditiile in 
care apar semnale perturbatoare cunoscute ori necunoscute si zgomote, sa faca totul 
cu comenzi mici (pentru a nu satura elementul de actionare), si totul in prezenta 
incertitudinilor de model... 
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Multe intr-adevar !!!! Este posibil sa facem toate aceste lucruri ? Daca DA, cum ? 
Incercam sa raspundem treptat acestor intrebari dezvoltand Teoria Sintezei. 


ATENTIE INSA: Modelele cu care lucram sunt black-box si nu putem umbla in interiorul 
lor (si oricum este nevoie de specialisti din domenii din cele mai diverse pentru a stii 
cum)! 


Cum procedam? Tratam modelul sistemului nominal ca atare (black-box) si anexam 
un alt sistem (numit uzual regulator si apartinand aceleiasi clase de modele) pe 
care il proiectam noi astfel incat sistemul rezultant format din sistemul nominal 
si regulator sa faca ceea ce dorim. 


Cum il “anexam” ? In serie, in paralel, in bucla de reactie pozitiva/negativa? 


Depinde ce dorim dar minimal dorim ca indiferent de sistemul original (stabil sau 
nu) sistemul rezultant sa fie stabil. 


Conexiunea Serie 


Fie P(s) = -+ pe care il inseriem cu un sistem proiectat de noi C' astfel incat sistemul 


~ 


rezultant sa fie stabil. Cum procedam ? Trebuie ca C sa fie de forma (s — 1)C (s) unde 
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C(s) sa fie strict stabila. Este sistemul rezultant strict stabil ? Da, dpdv intrare=iesire 
dar daca injectam un semnal intre cele doua sisteme NU va fi ! Mai mult, nu putem 
asigura un zerou EXACT in s = 1 si nici nu stim EXACT ca sistemul original are un pol 
in s=1 |! La cea mai mica abatere rezulta un sistem INSTABIL ! Deci nici un regulator 
nu-mi asigura in mod satisfacator cerinta de stabilitate (interna) ! 


Conexiunea Paralel 


Similar pentru conexiunea paralel ! Exemplificati voi cum ar trebui sa arate un regulator 
si de ce nici un regulator nu ne satisface | 


Cum mai putem lega doua sisteme in mod elementar ? 
Bucla de reactie (Feedback) 


Asa cum vom vedea printr-o analiza detaliata, conexiunea in bucla de reactie este 
CEA CARE REZOLVA CHESTIUNEA STABILITATII (interne) si, asa cum vom 
vedea, inca multe alte cerinte | 
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2. Bucla de Reactie (Feedback) 


Cel mai simplu sistem REAL cu bucla de reactie are trei componente (vezi Fig. 1): 


e Un sistem nominal ce este controlat (automatizat); 
e Un senzor ce masoara iesirea sistemului nominal; 


e Un sistem numit regulator (controller) ce genereaza intrarea sistemului nominal; 
(elementele de actionare sunt uzual incluse in sistemul nominal); 


Fiecare dintre cele trei componente ale buclei de reactie are doua 
semnale de intrare (unul din interiorul sistemului si unul  exogen) 
si un semnal de iesire. Semnalele au urmatoarele semnificatii: 

r: referinta (semnal de intrare) v: iesirea senzorului 

u: comanda (semnal de actionare) d: perturbatia externa 

y: iesirea sistemului (semnal masurat) n: zgomot de masura 
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Bucla de Reactie (Feedback) 


Cele trei semnale externe r, d si n se numesc intrari exogene. 


u Sistem 
Regulator narmina 
y Te 
| n 


Fig. 1: Sistem de control automat (elementar) 


Ce dorim ? y trebuie sa aproximeze o functie de r prespecificata si trebuie sa faca acest 
lucru in prezenta perturbatiei externe d, a zgomotului n si a incertitudinilor diverse de 
modelare (a sistemului nominal)! Mai mult, uzual nu acceptam nici comenzi u prea mari 
iar y este disponibil doar prin intermediul masuratorii v ! 


ipoteze matematice fundamentale (specifice TS 1): Semnalele din Fig. 1 sunt 
scalare, descrise de functii regulate sau singulare (generalizate) ce au transformate 
Laplace, cele trei componente ale buclei de reactie (sistemul nominal, regulatorul si 
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senzorul) sunt sisteme liniare invariante in timp avand functie de transfer rationala 
proprie (Totul este “civilizat” )! Deci “avem” frecventa si analiza urmatoare se face la 
nivelul transformatelor Laplace. 


Fiecare componenta din Figura 1 are doua intrari si o iesire putand fi exprimata sub 
forma (de exemplu in cazul sistemului nominal) 


„Pi |=[P P]| 4 | = Pid+ Pra 


Ipoteza aditionala: lesirile celor trei componente sunt functii liniare de suma (sau 
diferenta) intrarilor, i.e., ecuatiile sistemului nominal, senzorului si regulatorului sunt de 
forma 


y = P(d+u), 
v = F({y+n), 
u = Clr-—vw). 
Semnul “ — ” din ultima ecuatie este o chestiune de traditie (reactie negativa in 


(44 


electronica) putand fi inlocuit alternativ cu “+”. Diagrama bloc a acestor ecuatii este 
data in Figura 2 in care s-au omis semnele “+” in jonctiuni. 
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Fig. 2: Bucla de reactie (standard) 


Cand este (bine) definit asa ceva ? Toate functiile de transfer in bucla inchisa 
trebuie sa existe, i.e., functiile de transfer de la cele trei semnale exogene r, d si n la 
toate semnalele interne u, y, v si £1, £2 Si £3 | 


Prima concluzie: Este suficient sa analizam functiile de transfer de la r, d si n la z4, 
To Si £3 (celelalte rezulta din acestea). Scriind ecuatiile jonctiunilor obtinem 


ti = r-—/PHaz, 
t2 = d +CT, 
La = n+ Prs, 
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sau in forma matriciala 


1 0 F Tı r 
—C 1 0 Lo = d 
0 =P 1 La n 


Sistemul in bucla inchisa este deci bine definit daca matricea de dimensiune 3 x 3 
de mai sus este nesingulara, i.e., determinantul ei 1 + PCF nu este identic nul, i.e., 
1 +PCF ZO! In acest caz, functiile de transfer in bucla inchisa se obtin din 


=] 


Tı ji 0 F f 
T2 = -C 1 0 d 
Ta 0 -=P 1 n 
care este echivalenta cu 
Tı rl 1 -PF -F r 
te | =n C | —C F d l: (38) 
1 + PCF PC P 1 i 
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Problema: Am presupus intitial ca P, C, F sunt proprii (avand motivatie matematica 
serioasa). Rezulta si ca functiile de transfer in bucla inchisa (38) sunt proprii atunci cand 
sistemul in bucla inchisa este bine definit ? Din pacate NU fara ipoteze suplimentare ... 
III Introducem: 


Sistemul in bucla inchisa este bine definit in sens strict daca toate functiile de transfer 
in bucla inchisa sunt proprii (> 1+ PCF nu este strict propriu > 1 + PCF(œ) #0). 


Nota: Buna definire in sens strict implica automat si buna definire ! Este realist sa 
cerem unui sistem de reglare buna definire in sens strict ? DA !!! 


Argumente pro (si contra !): Functiile de transfer ale oricarui sistem sunt de fapt 
strict proprii (o sinusoida de amplitudine constanta si frecventa din ce in ce mai mare 
aplicata unui sistem conduce la o iesire ce tinde la zero). Da, dar la frecvente mari 
sistemele sunt in general foarte complexe si inceteaza in fapt sa fie liniare pe masura ce 
frecventa creste | De fapt avem nevoie ca unul singur dintre P, C sau F sa fie strict 
propriu si daca acest lucru nu se intampla automat la modelele cu care lucram oricum 
cerem |PCF(œ)| < 1 (motive profunde legate de stabilitate interna si “small-gain” ) 
ceea ce asigura automat buna definire in sens strict. 
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3. Stabilitatea Interna a Buclei de Reactie 


Problema (fundamentala): Ce notiune de stabilitate este adecvata sistemului din 
Figura 2 (stabilitate externa stricta, BIBO, asimptotica, Lyapunov etc.)? 


Raspuns: Niciuna nu este adecvata sistemului in bucla inchisa... 


Ce functii de transfer ne intereseaza de fapt ? In primul rand transferul |/O al sistemului 
in bucla r — y ! Dar chiar daca acest transfer este strict stabil este posibil ca alt transfer 
de la un exogen la un semnal intern al buclei sa nu fie stabil rezultand in distrugerea 
unei parti sau a intregului sistem !!! Deci avem nevoie ca toate functiile de transfer 


r, d, n —> Ti, T2, 3, V, Y, U, 


sa fie SIMULTAN stabile pentru a fi siguri ca totul este OK ! 


Exemplul 1. In Figura 2 fie C(s) = a, P(s) = a F(s) = 1. Functia de transfer 


r — y este strict stabila (in sens BIBO) dar functia de transfer de la d la y nu este. 
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Exemplul 2. In Figura 2 fie C(s) = a, P(s) = za F(s) = 1. Toate functiile de 
transfer exogen — semnal intern al buclei sunt strict stabile (in sens BIBO). 

Definitia 3. Sistemul in bucla de reactie din Figura 2 se numeste intern stabil daca toate 
transferurile in bucla inchisa sunt strict stabile (i.e., au polii in C_). 


Aparent trebuie sa verificam 18 functii de transfer dar de fapt sunt suficiente doar 9. 
Aratati de exemplu ca urmatoarele functii de transfer sunt suficiente pentru verificarea 
stabilitatii interne 

rd, n — 1, Xa, Ea. 


Este relativ complicat sa examinam chiar si 9 numitori si sa le calculam zerourile. Scriem 
fctiile de transfer rationale sub forma unor cături de polinoame coprime 


P= Œ, Ca, F= (39) 


si definim polinomul caracteristic al sistemului in bucla inchisa (din Figura 2) ca 
produsul numaratorilor + produsul numitorilor celor trei functii de transfer, 


x = NpNoNr + MeMcMr. 
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Poli: sistemului in bucla inchisa sunt zerourile polinomului caracteristic x. 
Teorema 4. Sistemul in bucla de reactie este intern stabil daca si numai daca polii 
sistemului in bucla inchisa sunt in C-. 


Demonstraţie. Pentru simplitate presupunem F = 1 (argumentatia in cazul general este 
similara dar este mai greu de urmarit). Din (38) avem 


Ti 1 1 -P -I T 
zə | = ———- C 1 —C d 
F3 L+ PC PG P 1 n 


si inlocuind (39) obtinem 


Tı 1 MpMc —NpMo —MpMo r 
mo | 2 l MN Mee -MeNe d (40) 
în NpNo+MpMo | NoNo NpMo  MpMo i 


in care polinomul de la numitor este exact polinomul caracteristic x. 
Suficienta: Sistemul in bucla inchisa este intern stabil daca polinomul caracteristic are 
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toate zerourile in C_. 

Necesitatea: Implica un argument mai subtil. Presupunem ca sistemul in bucla inchisa 
este intern stabil. Atunci toate cele noua functii de transfer in (40) sunt strict stabile (in 
sens BIBO), i.e., au poli numai in C_. Aceasta nu implica automat ca NpNo + MpMo 
are zerourile in C_ intrucat este posibil sa aibe un zerou in afara care sa fie zerou si al 
tuturor numaratorilor din (40) si deci sa se simplifice. 


In orice caz, polinomul x nu are un zerou comun cu toti numaratorii din (40) asa cum 
demonstram prin reducere la absurd. Fie so a.i. (NpNc + MpMo)(so) = 0. P.p. ca 
este zerou si al tutoror numaratorilor, in particular MpMo(so) = 0. Apar doua situatii 
posibile: i)Mp(so) = 0; ii) Mo(so) = 0. 

Cazul i) Mp(so) = 0; Automat Np(so) Æ 0 (Mp este coprim cu Np) si rezulta din 
elementul (1,2) al matricii ca Mce(so) = 0 si din elementul (3,1) ca Nc(so) = 0 ceea ce 
nu este posibil din nou din coprimitatea lui Na si Mo. 

Cazul ii) Mc(so) = 0. Automat atunci No(so) Æ 0 si rationamentul se repeta ca 
mai sus folosind elementele (2,3) si (3,1) ale matricii rezultand ca Mp(so) = O si 
Np(s9) = 0, ceea ce contrazice coprimitatea lui Mp si Np. 


Teorema de mai sus ne da o modalitate sintetica de verificat stabilitatea interna a unei 
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bucle de reactie: se verifica locatia in C_ a zerourilor unui polinom x (de exemplu 
calculand efectiv zerourile prin calculul valorilor proprii ale matricii companion asociate 
sau, mai simplu, folosind criteriul lui Hurwitz). 

Teorema 5 (Test de Stabilitate Interna). Sistemul in bucla inchisa este intern stabil 
daca si numai daca urmatoarele doua conditii sunt simultan indeplinite: 


a) Functia de transfer 1 + PCF are zerourile in C_; 
b) Nu au loc simplificari in Cu U Co cand formam produsul PCF. 


Demonstraţie. Sistemul in bucla inchisa este intern stabil daca si numai daca toate cele 


9 functii de transfer 
| FPF ep 


C 1 =C" 
Ie F 1 


1 
1+PCP 


sunt stabile. 

Necesitatea: Presupunem ca sistemul in bucla inchisa este intern stabil. Rezulta 
a a 1 3 a a a 

in particular ca por este stabila, deci 1 + PCF are zerourile in C_ ceea ce 


i | N N N 
demonstreaza necesitatea lui a). Fie P = sp, C = wo F= Me Din Teorema 4 
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rezulta ca polinomul caracteristic x = NpNcNr + MpMoMp are zerourile in C_. Prin 

urmare perechea (Np, Mc) nu are zerouri comune in Res > O si similar pentru toate 
. f NoNaN î 

celelalte perechi corespunzand lui PCF = pap. Deci si b) este adevarata ceea ce 

. | | P+HICIWp 

incheie necesitatea. 

Suficienta: Presupunem adevarate a) si b). Factorizam P,C, P ca mai sus si fie sg un 

zerou al polinomului caracteristic, i.e., 


(NpNeNp E îi MpMcMr)(s0) si (41) 


Aratam prin reducere la absurd ca Res < 0. Presupunem contrariul, Res > 0. 
Daca MpMcMpr(so) = 0 rezulta MpMoMr(so) = 0 ceea ce violeaza b). Deci 
MpMceMr(so) Æ 0 si impartind cu el in (41) obtinem 1 + MEME (50) — 0 de unde 
(1 + PCF)(so) = 0 ceea ce violeaza a). 


Concluziile “la zi” 


e Un sistem se comporta “civilizat” doar daca avem stabilitate interna, adica 
indiferent unde se injecteaza un semnal exogen de energie finita toate semnalele din 
sistem vor fi de energie finita; 
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e Stabilitatea interna nu se poate asigura cu conexiunile serie sau paralel ci doar 
cu bucla de reactie; 

e Stabilitatea interna se verifica pe baza faptului ca cele 18 functii de transfer din 
bucla de reactie sunt simultan strict stabile dar exista teste mult mai simple (vezi 
Teoremele 4 si 5); 

e Ramane sa investigam problema de sinteza: Dandu-se un sistem nominal oarecare P 
exista intotdeauna un regulator C care asigura stabilitatea interna ? Daca DA se cere 
constructia clasei regulatoarelor C' ce asigura stabilitatea interna. 
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4. Stabilizarea 


Consideram sistemul (mai simplu) cu bucla de reactie unitara din Figura 3. Aceasta 
diagrama arata situatia teoretica tipica in controlul automat (F = 1) in care se da 
P si se cere C ce satisface anumite cerinte: stabilizeaza intern, asigura performante, 
tolereaza incertitudini de modelare, limiteza comenzi, etc. Cu toate ca in general F Æ 1, 
dezvoltarile teoretice sunt destul de complexe chiar in acest caz mai simplu pentru a-l 
analiza separat. 

Ipoteza: P este strict propriu (bucla este automat bine definita in sens strict). 


Fig. 3: Sistem cu bucla de reactie unitara 
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Problema fundamentala: Dandu-se P exista intotdeauna un C care stabilizeaza intern 
sistemul in bucla de reactie? Daca DA, cati exista ? Cum ii putem caracteriza? 
Raspuns: DA, exista o infinitate de C care stabilizeaza intern si ei se pot da 
explicit prin intermediul parametrizarii lui Youla ! 


Parametrizarea lui Youla: Cazul P Stabil 


Deoarece cazul general este relativ complicat, incepem sa dezvoltam mecanismul necesar 
intr-un caz mai simplu: P stabil (P € S, unde © este inelul comutativ al functiior de 
transfer rationale proprii si stabile in sens strict). 
Teorema 6. Presupunem ca P € S. Clasa tuturor C pentru care sistemul in bucla de 
reactie din Fig. 3 este intern stabil este data de 


Tn QES}. (42) 


Teorema afirma doua lucruri: oricare ar fi C dat de (42) (cu un Q € S arbitrar) el 
stabilizeaza intern si, pe de-alta parte, pentru oricare Č care stabilizeaza intern se gaseste 
un Q a.i. O= o cu QES. 
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Demonstraţie. Presupunem ca Č stabilizeaza intern. Fie Q functia de transfer de la r 
la u 


C 
> 14, PU 
Evident Q e S si C = 5. 
Reciproc, fie Q e S. Definim O 
M 4 
C 1- PQ (43) 
Sistemul in bucla inchisa este intern stabil daca si numai daca cele 9 functii de transfer 
1 1 =P -1 
n C 1 -C 


sunt proprii si stabile. Substituind expresia (43) pentru C' aceasta matrice devine 


1- PQ -P(1-PQ) -(1- PQ) 
Q 1- PQ A fa (44) 
PO PO-PQ)  1-PQ 
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Cum Q € S si P € S este clar ca toate elementele acestei matrici sunt proprii si strict 
stabile. 


Observatia 7. Toate functiile de transfer in bucla inchisa sunt afine in parametrul liber 
Q, i.e., sunt de forma Ti + ToQ pentru anumiti Ti € S, 1» € S. Acest fapt simplifica 
major orice calcule necesare pentru alegerea parametrului Q a.i. sistemul in bucla inchisa 
sa satisfaca cerinte suplimentare (se poate interpreta ca o liniarizare). 


Trecem in continuare la analiza cazului general (P nu este neaparat stabil). In acest 
scop introducem factorizarile coprime peste S. 
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Factorizari Coprime 


Presupunem ca P nu este stabil si vrem sa gasim in continuare un C care stabilizeaza 
intern. Cum procedam ? 


Ideea de baza: Factorizam P peste S! 


De ce nu peste polinoame ? Sa presupunem ca scriem P = 3%, unde (N, M) sunt 
polinoame coprime. Atunci exista doua alte polinoame X si Y (ce se pot construi de 


exemplu pe baza algoritmului lui Euclid) ce satisfac identitatea lui Bezout 
NX +MY =1 


(ideea calauzitoare este sa folosim Teorema 1 care afirma ca sistemul in bucla inchisa 
este intern stabil daca si numai daca polinomul caracteristic nu are zerouri in Res > 0). 
Daca am seta pe baza identitiatii Bezout Č = — atunci polinomul caracteristic verifica 
NX +MY = l si deci nu are zerouri in Res > 0. De aici “concluzia” ca C stabilizeaza 
intern. 


Ce este gresit ? Y ar putea fi nul (deci C nu se poate construi astfel) sau OY < 0X 
(rezulta un regulator impropriu )! 
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Exemplul 8. Fie P(s) = 1 pentru care avem N(s) = 1, M(s) = s. O solutie a ecuatiei 
Bezout NX + MY = 1 este X(s) = 1 si Y (s) = 0 pentru care C := Š nu este definit. 
Alta solutie este X(s) = —s + 1, Y(s) = 1, pentru care C := X — —s + 1 nu este 


` Y 
propriu. 


Pentru a remedia aceasta situatie cerem ca M, N, X,Y sa fie elemente ale lui S si deci 
procedam la factorizari coprime peste S | 


Doua functii in $ se numesc coprime daca exista alte doua functii X,Y ambele in S 
a.i. sa aibe loc identitatea lui Bezout 


NAFTY E=] 


In particular, identitatea de mai sus implica ca N si M nu au zerouri comune in Res > 0 
sau la oo. Daca ar exista un astfel de sọ atunci am avea 


g= N (so) X (so) + M (so)Y (so) Z ls 


Conditia este si suficienta pentru coprimitate: daca N si M nu au zerouri comune in 
Re s > 0 sau la œ atunci sunt coprime (Demonstratia: Exercitiu ! ) 
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Factorizare coprima: Fie G matrice real rationala. O reprezentare de forma 
G = — N. MES, 


se numeste factorizare coprima a lui G peste S. Prezentam in continuare un mecanism 
(algoritm) de constructie a celor 4 functii din S ce satisfac cele 2 ecuatii: 


G = — NA FMY = 


Constructia lui M si N: Fie G rationala proprie cu coeficienti reali si fie G = o 


factorizare coprima peste polinoame (se obtine direct prin simplificarea oricaror factori 
polinomiali comuni). Evident On < Om =: k. Atunci G = & cu 


(4) m 


Ca (ul i 


este o factorizare coprima peste S a lui G. 
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Exemplul 9. Fie G = TÈ Atunci k := 2, n(s) = 1, m(s) = (s — 1)? si obtinem 


1)2' 
factorizarea coprima 
n(s) 1 m(s) _ (s— 1)2 


NS) = GI T GP MD 6+? 


Observatia 10. Asa cum se observa si in exemplul de mai sus, k este fixat de cerintele 
de coprimitate ale celor doi factori. Intr-adevar, un k mai mic creaza un M impropriu 
pe cand un k mai mare creaza doi factori M si N avand zerouri comune la infinit. 


Deci N si M se obtin facil. Obtinerea in continuare a lui X si Y este considerabil mai 
dificila necesitand algoritmul lui Euclid. 


Algoritmul Euclid calculeaza cel mai mare divizior comun a doua polinoame date n(A) 
si m(A). Cand n(A) si m(A) sunt coprime, algoritmul lui Euclid se poate folosi pentru 
calculul polinoamelor z(A) si (A) ce satisfac identitatea lui Bezout 


n(A)z(A) + m(A)y(Â) = 1. 
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Algoritmul lui Euclid 


Intrari: Polinoamele n(A) si m(A). Daca Om > On se interschimba rolurile lui m si n. 
Pasul 1: Imparte n la m pentru a obtine 


n = mqı + ru, Ori < Om 


(gradul restului mai mic strict decat gradul catului). 
Pasul 2: Imparte m la rı pentru a obtine 


Mm = Tri + To, rə < Orj. 


Pasul k: Imparte rg—ə la rgķ—ı pentru a obtine 
e-o = iki + Th Or < Orga. 


Opreste-te la Pasul k pentru care rą este o constanta nenula. 
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Ce am obtinut ? Scriind ecuatiile obtinute in forma compacta obtinem 


1 0 0 ri l —q 

d 1 0 To 0 1 N 

aar fefe e fla] 
—1 dk 1 Tk 0 0 


Prima matrice din membrul stang este inversabila si are inversa tot o matrice polinomiala 
(determinantul ei este o constanta) obtinandu-se resturile sub forma 


—1 


TI 1 0 0 1 —qı 

ra d 1 0 0 1 P 
aja na f fo ojia] 
Tk —] dk 1 0 0 


Deci r(A) := q(A)n(A) + q(A)m(A) in care q si q sunt polinoame functie de q;(A). Cum 
ry = cst. £ 0, fie z(A) := 1, y(A) := Hà care sunt polinoamele cerute. 
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Obtinerea Factorizarii Coprime peste S 


Am vazut algoritmul Euclid si ca rezultat identitatea Bezout pentru polinoame. Avem 
insa nevoie de identitatea Bezout peste functii din S (nu peste polinoame). 

Cum procedam ? Facem o transformare de variabila a.i. polinoamele in A sa genereze 
functii in S. 

Intrari: G 

Pasul 1: Daca G este stabila, setam N = G, M = 1, X =0, Y = 1 si STOP. Altfel 
continua; _ _ 

Pasul 2: Transforma G(s) in G(A) luand s := 134. Scrie G ca un cat de doua 
polinoame coprime 


GO) = 0) 


Pasul 3: Foloseste algoritmul Euclid (peste polinoame) si gaseste z(A) si y(A) a.i. 
nre d 0 NI 


Treci de la n(A), m(A), z(A), y(A) la N(s), M(s), X (s), Y (s) facand tranformarea inversa 
de variabila A = ==. 
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Parametrizarea lui Youla: Cazul General (P oarecare) 


Cu aceste preliminarii asigurate putem trece la constructia parametrizarii Youla pentru 
cazul in care P este arbitrar (dar propriu). 


Fie P = -> o factorizare coprima peste © si fie X si Y doua functii in S satisfacand 
identitatea lui Bezout 


NX +MY =1 (45) 


(toate acestea sunt furnizate de procedura elaborata anterior). 
Teorema 11 (Parametrizarea lui Youla). Multimea tuturor C pentru care sistemul in 
bucla de reactie din Figura 3 este intern stabil este data de 


X+MQ 
(FERo 


QES}. (46) 


Exercitiu: Aratati ca Teorema 6 se obtine prin particularizarea Teoremei 11 atunci cand 
PES! 


Demonstratia Teoremei 2 se face pe baza unei leme preliminare. 
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Lema 12. Fie P = — si Č = ME factorizari coprime peste ©. Sistemul cu bucla de 
reactie este intern stabil daca si numai daca 


(NN.+ MMo) e. 


Demonstraţie. Demonstratia lemei fiind aproape identica cu cea a Teoremei 5 este 


lasata drept exercitiu | 


Demonstratia Teoremei 11. Fie Q e S si definim 


c.. XEMQ 
IS 


Pentru a arata ca sistemul in bucla inchisa este intern stabil definim 


No:=X+MQ,  Mo:=Y-NQ 
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si incercam sa aplicam Lema precedenta. Pentru aceasta trebuie sa aratam succesiv ca 
C = ia este o factorizare coprima si ca (NN. + MM.) |! € S. Intr-adevar, avem ca 
NX + MY = 1 de unde rezulta automat ca NNo + MM. = 1 ceea ce demonstreaza 
prima parte a leoremei. 


Reciproc, fie C un regulator ce stabilizeaza intern bucla de reactie. Trebuie sa aratam 
ca exista Q € S a.i. 


X +MQ 
= 47 
C= FNQ (47) 
Fie C = Aei o factorizare coprima peste S si definim V := (NNo + MMo) 1. Din 
Lema 1 rezulta ca V € S. Avem ca 
NcV 
C= 48 
MoV Gi 


si incercam sa alegem Q a.i. (47) sa coincida cu (48) in sensul ca NoV = X + MQ 
si MoV = Y — NQ. II definim pe Q din prima ecuatie rezultand ca Q = SV si 


M 
calculand 
N(NoV — X) E MY + NX — NNcV 


Y-NQ=Y-— 
Q M M 
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1—NNOV  MMoV 
M M i 


arata ca o verifica de fapt si pe a doua. A ramas sa aratam ca Q astfel definit apartine 
lui S. Am vazut ca Q verifica 


MoV =Y — NQ, NoV = X +MQ. 
Scazand a doua ecuatie multiplicata cu Y din prima multiplicata cu X obtinem 
MecV X — NoVY = XY — XY —- (NX + MY)Q € Q = Nc VY — Mco VX 


si ultimul membru drept apartine lui $ intrucat toate elementele apartin lui S. Deci 
QES. 

Observatia 13. Analog ca in cazul in care P era stabila, toate functiile de transfer in 
bucla inchisa sunt functii afine de Q daca parametrizam C ca in enuntul teoremei in 
functie de Q. De exemplu functia de transfer de la r si d la € devine 


l P 
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Acest fapt va fi esential in sectiunile urmatoare in care parametrul Q se va determina 
a.i. sa satisfaca cerinte suplimentare de reglare (urmarirea asimptotica a unui semnal 
referinta si rejectia unui semnal de tip perturbator). 
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5. Performante Asimptotice 


Ce am obtinut pana in prezent ? Pentru o bucla de reactie standard (F = 1) 
am reusit caracterizarea clasei tuturor compensatoarelor care asigura stabilitatea interna 
(cea mai larga clasa ce ne intreseaza in aplicatii practice |) 


Ce mai dorim de la un sistem de reglare automat ? Minimal dorim sa-i asiguram 
o comportare dorita la anumite clase de semnale date si in conditiile in care pot aparea 
semnale perturbatoare si/sau zgomote. 


Ideal: Dorim ca la semnale precizate de intrare u € M iesirii sistemului sa-i corespunda 
un semnal dorit de iesire y € Y. Mai precis, dorim ca fiecarui u € M sa-i asociem un 
y € Y deci dorim ca sistemul sa reprezinte o aplicatie din U in Y precizata. 


Este realist sa cerem asa ceva ? In general NU !!! De ce ? Sunt mai multe motive: 

e u si y sunt functii de timp (sau rationale in domeniul transformatelor Laplace) iar 
sistemul este liniar, invariant in timp, descris de o transformata Laplace rationala proprie 
si stabila (cel putin in bucla inchisa). Deci pentru o pereche (uo, yo) fixata exista cel 
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mult un sistem T (rezulta din yọ = Tuo). Daca dorim ca la alta intrare u; sa corespunda 
o alta iesire yı rezulta in general alt sistem T' ! Concluzie: se poate cere asa ceva cel 
mult pentru o pereche de semnale ! 

e Este posibil asa ceva macar si pentru O SINGURA pereche de semnale? Nici 
macar ! Raportul = poate rezulta intr-un sistem instabil si deci neinteresant pentru 
aplicatii (semnalele de intrare interesante u sunt in general de tip persistent ca de 
exemplu treapta, sinusoide, rampe, etc.) 

e Mai mult, din raportul P pot rezulta sisteme improprii ! 


Ce este de facut? Relaxam cerintele pastrand esenta exigentelor! In general nu 
ne intereseaza ca un semnal de iesire sa aibe exact forma dorita de noi ci mai curand ne 
intereseaza ca dupa trecerea regimului tranzitoriu semnalul sa aibe o comportare dorita 
rezultand deci anumite cerinte de performanta asimptotica ! 


Clasa de semnale de intrare: semnale persistente (pentru noi = semnale cu transfor- 


mata Laplace rationala strict proprie et) in care polinomul u,;(s) are toate radacinile 
J 


antistabile, i.e., in afara lui C_). 


De ce facem astfel de ipoteze asupra semnalului? 
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e Transformata Laplace rationala: Lucram cu sisteme dinamice descrise de fctii de 
transfer rationale si deci este normal ca semnalele sa fie tot rationale (solutii ale unor 
ecuatii diferentiale ordinare in conditii initiale arbitrare). 


e Strict proprii: Semnalul in domeniul timp dorim sa fie functie in sens regulat (nu 
semnal generalizat). 


e Radacini antistabile: Putem considera semnale generale dar cele interesante din punct 
de vedere al proprietatilor asimptotice sunt doar cele antistabile. De ce ? Raspunsul 
unui sistem liniar y = Tu la o intrare u = u} + u— (u, partea antistabila si u_ partea 
stabila) este 


> ep te ae Su Yas u 


Deoarece sistemul este neaparat stabil (in bucla inchisa), rezulta ca y_ are toate radacinile 
in C_ si deci lim; sooy_(t) = 0 si dpdv asimptotic conteaza numai componenta 
antistabila. 


lpoteze: Presupunem configuratia clasica de reglare din Figura 4 
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Fig. 4: Sistem cu bucla de reactie unitara (F = 1) 


in care semnalele externe r si d sunt presupuse de tip persistent. Semnalul r numit 
referinta reprezinta semnalul pe care dorim sa-l urmareasca iesirea sistemului y iar 
semnalul d este semnalul perturbator pe care il dorim rejectat la iesirea sistemului 
(omitem momentan o discutie separata asupra lui n deorece F = 1 si deci presupunem 
ca este disponibila iesirea y (prim masurare exacta); mai mult, n fiind in general zgomot 
alb rejectia sa se poate trata doar cu instrumente specifice statisticii matematice) ! 
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Problema fundamentala a reglarii: Pentru configuratia de reactie din Figura 4 sa 
se gaseasca un sistem (numit compensator regulator) a.i. sistemul in bucla inchisa sa 
indeplineasca simultan cerintele: 


Stabilizare (S): Sa fie intern stabil, adica C sa fie un compensator stabilizator; 


Reglare (R): Sa indeplineasca cerintele de reglare asimptotica adica 


lim e(t)=0 WER, 


t— 00 


lim e(t)=0 vdeP, 


t— oo 

adica C sa fie un compensator regulator pentru clasa referintelor R si clasa pertubatiilor 
D. 

Observatia 14. e Cerinta de reglare (performanta asimptotica) nu are sens fara 
asigurarea simultana a conditiei de stabilizare ! Nu are deci sens sa abordam separat 
problema reglarii de cea a stabilizarii de exemplu prin parametrizarea clasei tuturor 
regulatoarelor si ulterior intersectarea cu clasa celor care stabilizeaza. Problema se 
abordeaza prin deducerea subclasei compensatoarelor stabilizatoare ce asigura in 
plus si reglarea ! 
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e Cerinta de reglare expliciteaza de fapt doua cerinte individuale de comportare la doua 
clase de semnale de intrare: clasa de referinte si clasa de perturbatii sub forma: 
Urmarirea referintei (Rr): 


lim e(t) =0, YvreR ,„d=0 


t— 00 
Rejectia perturbatiei (Rd): 


lim e(t)=0, vdeP „r=0. 


t— 00 


e Ambele clase de semnale R,D sunt presupuse cunoscute (date prin cerinta de 
proiectare). Cele doua cerinte se considera separat rezultand in final si compunerea 
efetcelor lim; „x c(t) = 0,Yr e R,vd e D (acest fapt este o consecinta a principiului 
superpozitiei aplicabil sistemelor liniare). 

e Cerinta de urmarire a referintei se poate scrie echivalent 


Jim y(t) = lim r(t) YrER ,d=0. 


t— 00 
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Reformularea problemei reglarii 


Fie P = — si C = A factorizari coprime ale lui P si respectiv ale lui C' peste S si fie 


r = a, d = % factorizari coprime peste polinoame. Deorece r si d sunt presupuse 


Tm dm 


semnale persistente avem automat ca 
Ora < Olm Odla < Odi 


si radacinile lui rm si dm sunt in C, U Co. 
Expresia lui € ca functie de semnalele de intrare d si r este 


l P MMe Fa N Me (lg, 


E ea aa i i ll 49 
14 PC 1+ PO NN. + MMrm NN. + MMedm (49) 


€ 


Teorema valorii finale: Daca e(s) este o transformata Laplace rationala fara poli in 
Re s > 0 cu posibila exceptie a unui pol simplu in s = 0 atunci lim, sos €(t) exista si se 
poate calcula cu alternativ prin 


Jim e(t) = lim se(s). 
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Problema reglarii capata atunci urmatoarea formulare echivalenta: 


Problema reglarii: Sa se gaseasca doua elemente in S notate N. si Me, cu - proprie, 
C 
care satisfac simultan cerintele: 


Stabilizare (S’): NNO MTI CS 


Reglare (R'): Polii lui se(s) C C_ si 


lim se(s) = OflyreRr,a=0, lim se(s) = Olvaep,r=o- 
s—0 


s= 


Observatia 15. Chiar in aceasta formulare problema este relativ complicata intrucat cele 
doua conditii (S') si (R') trebuie asigurate simultan ! Asa cum vom vedea putin mai 
tarziu, problema se simplifica considerabil daca folosim pentru Ne si Me expresiile deduse 
la stabilizare, adica daca scriem formula parametrizata a compensatoarelor stabilizatoare 
si impunem sa satisfaca exclusiv conditia (R'), (S') fiind automat satisfacuta in acest 
caz. |n aceasta idee este formulata si teorema reglarii din sectiunea urmatoare. 
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0. Solutia Problemei Reglarii 


Teorema 16. Presupunand conditia (S) (sau echivalent (S')) din problema reglarii 
indeplinita, conditia (R) (sau echivalent (R')) este indeplinita daca si numai daca 
urmatoarele doua conditii sunt simultan indeplinite: 


1 1 
x — E 5 Z(rm) C ZI 


| 1 
OPE E TI PO) Zr) C P(P)U PU) 


) o Z(dm) C Z(P)U P(C) 


Demonstrație. Sa observam intai ca sirurile de echivalente din enunt rezulta automat 
din conditia ca rm Si dm au toate zerourile in C+ U Co si sistemul in bucla inchisa este 
intern stabil. 


Conditiile sunt suficiente: Intr-adevar, conditiile fiind indeplinite rezulta din (49) ca 
e(s) are toti polii in C_. Mai mult, deoarece r si d sunt descrise de rationale strict 
proprii rezulta si ca toti polii lui se(s) sunt in C_. 
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Deoarece S$ x — E€ Ssi PS x + € S avem ca s nu este pol al acestor rationale si din 
(49) rezulta 


lim se(s) = 0. 


Deci conditia (R') este in totalitate indeplinita. 


Conditiile sunt necesare: Deorece (R') este indeplinita, rezulta din (49) ca simultan 
i e Dl RE Uli (50) 
r m 


Cum sistemul este intern stabil si deci PG e S iar Z(rm) C C+ U Co rezulta automat 
din conditia (R') folosita in prima relatie (50) ca rm trebuie sa dispara complet prin 
simplificare. Rationamentul pentru a doua conditie din (50) este similar rezultand in 
final necesitatea conditiilor (i) si (ii) din enunt. 


Observatia 17. e Sa observam ca PG are un zerou instabil in sg (care sa simplifice 
corespunzator “modelul” semnalului referinta continut in rm ) daca si numai daca functia 
de transfer in bucla deschisa L = PC are un pol acolo. 

e [eorema de mai sus se mai numeste si principiul modelului intern: pentru a avea 
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urmarire asimptotica a semnalului referinta trebuie ca polii semnalului referinta sa apara 
ca poli in modelul intern al buclei deschise L. 
e O analiza asemanatoare se poate face si in privinta semnalului perturbator d. 


Exemple de aplicare: 


Sa presupunem ca avem un sistem P si dorim sa construim un regulator care sa urmaresca 
un semnal de tip treapta sau de tip rampa 


r(t) = c, daca t>0, fidea ct, daca t>0, 
"| 0. daa t<0, ? 10, daca t<0, 


unde c € R*. Alternativ, dorim sa rejectam de exemplu o perturbatie de tip sinusoidal, 


d(t) = sin(ut)1(8). 


Exemple tipice pentru referinta treapta il constituie reglarea automata a vitezei unui 
automobil (cruise control) sau reglarea temperaturii dintr-o camera cu ajutorul unui 
termostat ce comanda centrala termica. Un exemplu tipic pentru reglarea la rampa il 
constituie o farfurie de radar sau de TV prin satelit ce trebuie sa urmareasca orbita 
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circulara a unui satelit (negeostationar). Un satelit ce se misca pe o orbita circulara cu 
o viteza unghiulara constanta acopera un unghi ce este functie liniara de timp (adica 
o rampa). Perturbatii sinusoidale apar in medii magnetice sau bune conducatoare de 
electricitate acolo unde este folosita de exemplu tensiunea alternativa. Aplicand teorema 
reglarii in aceste cazuri obtinem urmatoarea lema. 

Lema 18. Prespunem ca sistemul in bucla inchisa este intern stabil si d=0. 

(a) Daca r este o treapta atunci are loc reglarea (i.e. lim; „os e(t) — 0) daca si numai 
daca S := pg are cel putin un zerou in origine (s = 0). 

(b) Daca r este o rampa atunci are loc reglarea (i.e. lim; „os e(t) — 0) daca si numai 
daca S := PG are cel putin doua zerouri in origine (s = 0). 

Prespunem ca sistemul in bucla inchisa este intern stabil si r = 0. 

(c) Daca d este o sinusoida d(t) = sin(wt)l(t) atunci are loc rejectia perturbatiei (i.e. 
lim sos €(t) — 0) daca si numai daca P are un zerou in s = jw sau C are un pol s = jw 
(automat va avea si s = — jw deoarece sistemele au coeficienti reali). 


Demonstraţie. Rezulta direct din aplicarea transformatei Laplace asupra lui r(t) si d(t) 
si Teorema reglarii. 


Exemplul 19. Sa luam un exemplu banal. Fie P(s) = +, C(s) = 1. Functia de transfer 
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de la r la € este 
S 


1l+s-l epil 

Prin urmare polul din bucla deschisa din s = 0 devine zerou al functiei de transfer de 
la referinta la eroare si acest zerou anuleaza polul lui r(s) = 1 rezultand un e(s) fara 
poli instabili. Deci acest sistem asigura in bucla inchisa urmarirea unei referinte de tip 


treapta unitara. 


Tey = 


Cu toate ca teorema reglarii ne da un instrument teoretic extrem de util in abordarea 
analitica a solutiei problemei reglarii, mai sunt inca de facut cativa mici pasi in explicitarea 
clasei compensatoarelor stabilizante si gasirea parametrului Q care satisface conditiile 
enuntate. 


Clasa compensatoarelor care asigura cerinta primordiala de stabilitate interna este data 


de C = a in care Q este arbitrar in S si P = — si NX + MY = 1, cu toate 


N.M,X,Y €S. 


Functiile de transfer ce intervin in Teorema reglarii sunt date de 


1 M Me N Me 


CD a Sa l Ta = — PO = -c 
Lr NNI+-MM ~ “d NN. + MM. 


CAPITOLUL 4: Concepte Fundamentale ale Buclei de Reactie 247 Solutia Problemei Reglarii 


care devin dupa inlocuirea formulelor pentru Me := Y — NQ si Ne := X + MQ, 


Te = M(Y —- NQ), Tea = N(Y — NQ). 


Scriind acum rm Si dm in forma 


Tm = k„(s = Srı)(s — Sr2) E (s == Srk) 


dm(s) = kals — sai) (s — 842)... (s — Sae) 
unde kr, ka € R si radacinile Sri, Sq E C+ U Co, t = 1,...k,j = 1, 


complex conjugate. 


..., é apar in perechi 


Cu aceste notatii, conditiile de interpolare enuntate in Teorema reglarii devin pentru 


urmarirea referintei 


M(3r)(% (sre) — N(8,4)Q(s74)) = 0, 
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(51) 


Solutia Problemei Reglarii 


iar cele pentru rejectia perturbatiei devin 


| | (52) 
N (sae) (Y (sae) — N (sae)Q(sae)) = 0. 


e |n ecuatiile de mai sus N, M, X,Y sunt elemente cunoscute (determinate anterior) in 
S si deci valorile lor in punctele sp; Si Sp; sunt numere reale (sau complexe) ce se pot 
calcula direct. Necunoscuta o constituie parametrul liber Q € S care se determina a.i. 
conditiile de interpolare de mai sus sa fie indeplinite. 

e Daca s,; sau sq; sunt reali, ecuatia corespunzatoare va avea toti coeficientii reali 
si evident si valoarea rezultata pentru Q(sr;) sau Q(sq;) va fi reala. Daca radacina 
corespunzatoare este complexa atunci automat valoarea pentru Q(sri) sau Q(sq;) va fi 
complexa rezultand doua ecuatii corespunzatoare partii reale si celei imaginare avand ca 
solutii Re Q(-) si Im Q(-). Ecuatiile care rezulta pentru radacina complex conjugata Sr; 
sau 5r; vor fi identice cu cele obtinute pentru Sp; sau Spj si nu mai este necesar sa le 
rezolvam/scriem. 

e Parametrul Q este complet arbitrar in © insa conditiile de interpolare impun esential- 
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mente conditii asupra zerourilor sale. In consecinta pentru simplificarea cat mai mult 
posibil a calculelor necesare se ia un Q € S de o forma cat mai simpla de tipul 


Q(5) = ap Fais+... tirgi ~ 
(s + 1) 
sau 1 
Q(s) = 40 + “a T I -+ sei + hi I a PRE 


Numarul minim de parametri cu care se incearca gasirea lui Q este dat de numarul de 
conditii de interpolare ce trebuie satisfacute, i.e., k + £ si fixeaza gradul (minim) al 
numaratorului si numitorului lui Q la &k+/— 1 (Q trebuie sa fie propriu). Se observa ca 
pentru orice valori ale parametrilor ao, a1,..., Grup 1 avem Q €E S. 

e Daca modelul intern este asigurat de sistemul original P atunci ecuatiile (51) 
si (52) sunt satisfacute indiferent de valoarea lui Q. Intr-adevar, daca in cazul 
urmaririi referintei avem ca s,. este pol al lui P, atunci automat M (si) = O si ecuatia 
respectiva 

Misa Y (Sy) — N (oi 0 as] ] = 0, 


nu conduce la nici o restrictie asupra lui Q(s,.;). Similar pentru cazul rejectiei perturbatiei 
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daca sq este in modelul original P(s) sub forma unui zerou atunci ecuatia 


IN (saj)(Y (sa) — N (saj)Q(saj)) = 0, 


este automat satisfacuta pentru ca N (saj) = 0 si deci nu apare nici o restrictie asupra lui 
Q)(s3qj). Aceste fapte ne permit reducerea complexitatii parametrului Q (si implicit 
a complexitatii regulatorului rezultant) in cazurile mentionate. 

e Daca apare necesitatea asigurarii unui model intern cu radacini multiple ecuatiile 
(51) si (52) trebuie completate cu variantele lor derivate de k — 1 ori, unde k este 
multiplicitatea radacinii lui r(s) sau d(s). 

e Ecuatiile (51) sunt satisfacute sau daca M (spi) = 0 caz in care ecuatia se elimina si 
se scade cu 1 sau cu 2 gradul lui Q (caz real sau complex) sau daca 


Dupa eliminarea tuturor ecuatiilor triviale (generate de prezenta chiar si partiala a 
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radacinilor in modelul intern) rezulta deci ca Q se alege a.i. sa satisfaca ecuatiile 


ceea ce revine la gasirea parametrilor ag, a1, ...a@k+e—1 ce satisfac sistemul de ecuatii 


Y (sri) _ 1 1 P 
a i ao F cd ES] EE . e n zi Ak+L—1 g 41) EFT? Vi — 1; e.. k, 
Sdj) _ 1 1 -o 
N (saj) = ao + alagi + e.. + CR Tsy IET Vi > o... EA 
Acest sistem este un sistem liniar de k +£ ecuatii si k +£ necunoscute ag, a1, ..., ak+£—1, 


avand matricea sistem de tip Vandermonde. Notam in continuare unitar cu s; punctele 
de interpolare atat corespunzatoare referintei cat si perturbatiei (mai exact s; := Sri ptr. 
i = 1,...,k Si 8; := Sasu ptr i = k+ 1,...,k+ 4. Sistemul Vandermonde se rescrie 
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Aa, = b unde 


1 dle e IE: ENI Y (s1) 

sl (s1+1)2 “ri(s N (s1) 

O a ind ici, AA ao Y (s2) 

(s2+1)? "(so I)kit-l a N(s2) 

.— i i .— l .— | 
A:= | | aa := b:= 

RE N EI NI Opel Yena 
Sktl (sktl)? 7 (sgpet1)EtE1 N(sk+e) 


Sistemul are solutie unica daca si numai daca matricea A este inversabila sau det A Æ 0. 
Se stie ca pentru o matrice Vandermonde are loc formula 


det A = liae — Sj) 


si deci problema are solutie unica daca si numai daca punctele de interpolare sunt 
distincte. 

e Daca punctele de interpolare nu sunt distincte (apar radacini s; cu ordinul de 
multiplicitate m > 1 atunci se scrie ecuatia de baza si inca m — 1 derivate ale sale. 
Matricea A devine de tip Vandermonde extinsa avand cele m linii corespunzatoare lui s; 
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de tipul 


| că 1 1 
ae aP E aa i 
0 1 TES eri (k +t- LTD 
0 0 2 na (+ 6 — a 
si se actualizeaza corespunzator coeficientul termenului liber He ce corespunde com- 


ponentelor derivate. 


e Din cele de mai sus aparent reiese ca problema reglarii are intotdeauna solutie. Acest 
lucru nu este adevarat, si rezulta din faptul ca este posibil ca N(s,) = 0 caz in care 
automat M(s;) Æ 0 iar ecuatia de interpolare M(s,)(Y(s;) — N(s;)Q(s;)) = 0 poate 
fi satisfacuta daca si numai daca Y(s,) = 0. Aceasta egalitate nu poate insa avea loc 
intrucat rezulta ca 
(NX + MY)(sn)=0 #1. 

Teorema urmatoare sintetizeaza conditiile necesare si suficiente de existenta a unei solutii 
a problemei reglarii. 

Teorema 20. Fie P = — o factorizare coprima peste S si X,Y e Sai. NX+MY =. 
Fie r = si d = go semnale persistente de tip referinta si respectiv de tip perturbator 
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factorizate coprim peste polinoame. Fie rı si dı c.m.m.d.c. (cel mai mare divizor comun 
— peste polinoame — ) al perechii (rm, M) si respectiv al perechii (dm, N) si fie 


PSD > Tm d, _ d Aa , si sad 
Tm = ; m = , 
TI di 


Problema reglarii are o solutie daca si numai daca N nu are nici un zerou comun cu rd. 


Demonstraţie. Conditia este necesara: P.p. prin absurd ca s; este un zerou comun 
atat al lui N(s) cat si 7d(s). Dupa simplificarea c.m.m.d.c. 7 si d, in ecuatiile de 
interpolare rezulta ca s; trebuie sa satisfaca ecuatia 


Y (s:) — N(s:)Q(s:) = 0 


care asa cum am vazut deja este posibila daca si numai daca N(s;) Æ 0. De unde 
contradictia. _ 

Conditia este suficienta: Intr-adevar, daca rd nu are zerouri comune cu N rezulta ca 
N (si) Z 0 pentru fiecare zerou al lui rd(s) si deci putem scrie un sistem Vandermonde 
corespunzator care are intotdeauna o solutie unica. 
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Procedura pentru rezolvarea problemei reglarii 


Pasul 1: Se rezolva problema stabilizarii folosind parametrizarea lui Youla pentru clasa 
tuturor compensatoarelor stabilizante; 

Pasul 2: Se analizeaza specificatiile de performanta asimptotice si se reduc la conditii 
de interpolare asupra parametrului Q; se trateaza distinct cazul cu radacini multiple si 
se analizeaza existenta solutiilor; 

Pasul 3: Daca problema are solutie se rezolva sistemul Vandermonde corespunzator si 
se determina parametrul Q care satisface deci conditiile de interpolare; 

Pasul 4: Se inlocuieste parametrul Q in expresia parametrizata a compensatorului si se 
obtine compensatorul regulator. 

Observatia 21. Toate dezvoltarile de mai sus se pot face folosind modele polinomiale 
insa natura intima a problemelor de reglare (cu asigurarea prealabila a stabilizarii) indica 
ca mult mai potrivita abordarea de mai sus ce foloseste factorizari coprime peste S | 
Observatia 22. Se poate formula o problema conexa problemei de reglare deja formulate, 
numita problema de reglare structural stabila. Aceasta cere ca pentru o vecinatatate 
intreaga (cat de mica) a modelului nominal compensatorul sa stabilizeze si sa continue 
sa regleze. Se poate arata relativ usor ca acesta problema are solutie daca si numai daca 
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cel mai mic multiplu comun al lui rm si dm notat rd divide Y — NQ adica daca si numai 


daca PI 
N(s,) 740, vs ai. rd(s;)=0 


& rd(s) nu are zerouri comune cu N(s). 


Acest fapt implica ca modelul semnalelor externe (exogene) r si d trebuie sa fie integral 
inclus in modelul regulatorului si nu poate fi folosita (in caz ca exista) partea comuna 
cu modelul nominal P(s). 


Din punct de vedere al aplicatiilor concrete ale reglarii, avand in vedere ca modelul 
nominal este aproape intotdeauna imprecis se impune deci includerea INTEGRALA a 
modelului exogenului in regulator (chiar daca este deja total/partial continut in P). 


In general aceasta masura de precautie nu este suficienta pentru multe aplicatii intrucat 
nu garanteaza rezolvarea problemei decat pentru o vecinatate arbitrar de mica a modelului 
nominal. Formularea corecta si rezolvarea completa a problemei se face in cadrul teoriei 
robustetii in care se gasesc acei () care maximizeaza vecintatea in care are loc stabilizarea 
si/sau performanta asimptotica. 
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Alte cerinte de performanta realizabile cu Q 


Parametrul Q poate fi folosit pentru realizarea altor cerinte de performanta: 


e Obtinerea anumitei comportari dinamice dorite a sistemului in bucla inchisa prin 
asigurarea anumitor locatii a polilor in bucla inchisa) — (exercitiu: propuneti o metoda 
pentru aceasta !); 


e Asigurarea stabilizarii robuste ce permite cele mai mari incertitudini (intr-o anumita 
norma) asupra modelului nominal P cu pastrarea stabilitatii interne; 


e Performanta robusta ce permite asigurarea proprietatilor de tip asimptotic pentru o 
incertitudine cat mai mare in modelul nominal P; 


e Asigurarea reglarii cu o clasa de regulatoare ce permit limitarea comenzii la o anumita 
valoare prescrisa; 


e Sa se stabilizeze un sistem stabil cu un compensator stabil; 
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e Problema stabilizarii simultane a doua sau mai multe sisteme nominale date folosind 
acelasi regulator 


e Combinatii ale tuturor sau doar unora dintre cerintele de mai sus. 


Desigur, combinand toate sau chiar o parte dintre aceste cerinte problemele se complica 
major si exista desigur situatii in care nu se pot asigura simultan toate cerintele. Un studiu 
mai detaliat al conditiilor in care anumite cerinte din cele de mai sus pot fi satisfacute se 
va face in semestrul al II-lea (implica introducerea unui formalism matematic adecvat). 
Desigur exista si numeroase probleme deschise (a caror solutie nu a fost inca data) pe 
care sunteti invitati sa le incercati ! 
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Exemple de Probleme (tipice de examen) 


Exemplul 23. Fie un sistem nominal stabil P. Sa se gasesca toate regulatoarele care 
urmaresc o referinta treapta. Mai mult, sa se precizeze subclasa solutiilor structural 
stabile. 

RRN Orice reglare se face cu stabilizare prealabila. P fiind stabil, rezulta ca 
C= 750 da clasa tuturor compensatoarelor stabilizatoare cand Q parcurge S. In 
cazul in care P este stabil, avem trivial N = P, M = 1, Y = 1, X = 0. Pentru 
un semnal de tip treapta r(t) = 1(t) avem r(s) = 1 si deci rm(s) = s cu singurul 
zerou in s = 0. Cum M(0) = 1 Æ 0 rezulta ca singura varianta de a satisface ecuatia 
M(0)(Y (0) — N(0)Q(0)) = 0 este sa avem P(0) = N(0) 740. In concluzie problema 
are solutie daca si numai daca s = 0 nu este zerou al lui P caz in care orice Q € S cu 


Qo = B(aj satisface cerinta de reglare, i.e., 


Aceasta expresie da si clasa solutiilor structural stabile (explicati de ce !!!). 
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Exemplul 24 (Exemplu detaliat la seminar). Fie P(s) = EE Sa se gaseasca 
un regulator ce asigura urmarirea unei referinta treapta r(t) = 1(t) si rejectia unei 
perturbatii sinusoidale d(t) = sin10t. Rezolva orice regulator gasit problema structural 


stabila ? 


Schita de rezolvare: Trebuie gasita intai clasa compensatoarelor stabilizante. P fiind 


X+MQ 


instabila aceasta are forma C = Şi 


= lR X 


PASER S] SEN 


z cu n(A) = 17, 


Aplicam algoritmul lui Euclid pentru n(A) si m(A) obtinand succesiv 


1 5 1 
q(A) = G rı(à) = z^ Tp 

36 1 6 
q2( A) = ZA 6’ Trə(A) = 25 
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cu Q arbitrar in S. Deci trebuie intai gasiti 
N,M,X,Y folosind procedura de la pagina 29. Mai intai calculam 


m(A) = 6åf — 5A + 1. 


Solutia Problemei Reglarii 


si ne oprim intrucat am obtinut un rest constant nenul. Ecuatiile verificate sunt 


N = mgi FT 
m = Tiqa + r2 


de unde 
rə = (1 + q1q2)m — qon. 
Deci luam i 
PRR PE pe e N oma a 
r2 r2 
Revenind in variabila s prin transformarea inversa À = = obtinem 
1 (s — 1) (s — 2) 
N = ~ M — 
(5) (s +1)? (4) (s + 1)? 
l 19s — 1 1 s+ 6 


Deci am deteminat toti C care asigura stabilitatea interna. 


CAPITOLUL 4: Concepte Fundamentale ale Buclei de Reactie 262 Solutia Problemei Reglarii 


a În ai __ 100 —. = 
Avem r(s) iri d(s) = -2100 de unde rm(s) = s si dml8) S 100. Observam ca 
nici o radacina a acestor doua polinoame nu apare ca pol respectiv ca zerou al sistemului 
initial P si deci problema de reglare structural stabila va avea aceeasi solutie cu cea 


obisnuita data de cele doua conditii de interpolare ce trebuie satisfacute: 


M(0)(Y (0) — N (0)Q(0)) = 0, 
N(105)(Y (103) — N (105)Q(105)) = 0. 


Inlocuind valorile corespunzatoare ale lui N, M, Y obtinem echivalent 


QOJ) = FAR = -94+ 10j & ReQ(10j) = —94, ImQ(10j) = 70. 


Avem de satisfacut trei conditii de interpolare si deci vom lua un polinom Q (in variabila 


+) cu trei coeficienti deci de gradul 2: 


1 1 
Q(s) =y AE -+ BEET 
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Ramane sa determinam coeficientii ag, a1, a2 pe baza conditiilor de interpolare. Re- 
zolvand sistemul liniar de tip Aa, = b ce se obtine rezulta coeficientii ap = —79, 
aı = —T23, a2 = 808. De aici rezulta expresiile 


__—7952 — 881s + 6 


Q(5) __—60s* — 598s? + 2515s? — 1794s + 1 
a a 


C(s) = s(s2 + 100)(3 +9) 


care reprezinta solutia cautata. 

Exemplul 25. Sa presupunem ca in exercitiul precedent se cere realizarea unei urmariri 
a unei referinte de tip rampa, in contextul in care perturbatiile nu apar (sunt nule). 
Semnalul rampa are transformata Laplace 


si deci rm(s) = s7. Conditia de interpolare se pune astfel incat M(s)(Y(s) — N(s)Q(s)) 
sa aibe un zerou dublu in s = 0. Cum M(0) Æ 0 rezulta echivalent ca Y(s) — N(s)Q(s) 
trebuie sa aibe un zerou dublu in s = 0 si deci Q(s) trebuie sa fie de gradul 1 pentru a 
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avea doi parametri liberi ag si aj, i.e., 


din care rezulta sistemul corespunzator cu solutia ap = 13,a4 = —7. Similar se scrie 
acum Q si respectiv C. 
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